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SUMMARY 

After  a  brief  presentation  of  the  salient  observational  facts  about 
tropical  cyclones,  the  two  principal  models  for  the  cyclone-scale  structure 
and  maintenance  of  the  quasi  steady  mature  stage  are  described-- the  older 
model  by  Kiehl  and  Malkus  and  the  recently  developed  model  by  Carrier.  One 
point  of  disagreement  between  the  two  models  is  that  Riehl  and  Malkus 
postulate  that  greatly  augmented  sea-to-air  transfer  of  latent  and  sensible 
heat  is  required  to  sustain  a  hurricane,  while  Carrier  states  that  (because 
the  storm  ccnvects  a  large  mass  of  warm  moist  air  with  it)  the  local  sea-to- 
air  transfer  of  total  enthalpy  is  about  the  same  whether  a  hurricane  is 
present  or  not.  Analyses  relating  to  maximum  swirl  speed  estimates  and  to 
surface-frictional-layer  dynamics  and  energetics  are  reviewed  to  support 
the  Carrier  model.  Reservations  concerning  both  the  frictional  inflow  layer 
modeling  and  also  the  interpretation  of  temperature  measurements  by  Riehl 
and  Malkus  are  stated  in  detail.  Since  the  cyclone-scale  structure  of  a 
**  fully  developed  tropical  cyclone  is  now  believed  to  be  correctly  sorted 

out  by  Carrier's  model,  attention  is  then  turned  to  refinements  of  Carrier's 
earlier  outline  of  tropical -cyclone  intensification,  during  transition 
from  depression  to  hurricane.  Intensification  theory  necessitates  careful 
consideration  of  the  problem  that  has  preoccupied  most  numerical  modelers-- 
parameteri zation  of  cumulus  convection.  Here  some  preliminary  and  tentative 
modeling  of  the  development  of  a  warm  core  is  attempted,  in  which  organized 
convection  lifts  and  ejects  initially  present  air,  whose  thermodynamic  state 
is  maintained  by  relatively  slow  ambient  processes  (cumulus  convection, 
turbulent  mixing,  radiational  cooling).  More  basic  incorporation  of  Charney's 
concepts  concerning  "conditional  instability  of  the  second  kind"  on  the  cumulus 
scale,  just  as  Carrier's  model  in  a  sense  already  utilizes  the  concepts  on  a 
cyclone  scale,  is  cited  as  probably  the  source  of  the  next  important 
progress  in  modeling  of  tropical  cyclone  intensification. 
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I.  INTRODUCTION 


I. A  MOTIVATIONS  FOR  HURRICANE  MODELING 

In  an  average  year  the  Atlantic  and  Gulf  coast  states  suffer  over 
$100  million  damage  and  50-100  fatalities  owing  to  hurricanes;  in  a  severe 
year  damage  will  exceed  $1  billion  (Meyer  1971).  Hurricane  Camille  alone 
in  August  1969  cost  about  500  lives  and  $1  billion  in  destruction 
(De  Angel  Is  1969).  The  threat  comes  in  the  form  of  wind  [over  200  mph  is 
known  (Lear  1969)];  in  rainfall  [27  In.  in  24  hours  (Schwarz  1970)-since 
1886  hurricanes  have  caused  over  60  floods  in  the  U.S.  (Alaka  1968)];  and 
in  storm  surges  [coastal  ocean  levels  have  risen  as  much  as  15-20  feet 
(Alaka  1968)].  On  the  positive  side  Is  the  fact  that  on  a  long-term  basis 
the  average  hurricane-associated  rainfall  over  the  Eastern  states  is  a 
substantial  fraction  (perhaps  onr:-th1rd)  of  the  total  rainfall;  without  it, 
droughts  seem  Inevitable. 

The  threat  to  execution  of  military  missions  Is  evident;  there  are 
annually  about  fifty  hurricanes  on  a  global  basis,  affecting  all  oceans 
except  the  S.  Atlantic  (Palmdn  and  Newton,  1969).  Actually  unnecessary 
preparations  by  DOD  Installations  owing  to  false  hurricane  and  typhoon 
warnings  is  annually  about  $8.3  mllllon-this  Is  aside  from  diverslon-of- 
manpower  costs  (Malone  and  Lelmer  1971).  In  path  predictions  for  either 
military  or  civilian  use,  three  times  the  area  actually  hit  by  a  tropical 

cyclone  is  typically  placed  under  hurricane  warnings  (Meyer  1971;  Malone 
and  Lelmer  1971). 

These  considerations  alone  would  justify  an  effort  at  understanding 
the  genesis,  steady-state  structure,  and  decay  of  tropical  cyclones  to  be 
able  to  forecast  their  onset,  Intensity,  and  path.  However,  there  Is  also 
an  additional  motivation  of  partlci  lar  relevance  to  the  current  study. 

A  hurricane  has  been  estimated  to  have  a  kinetic  energy  comparable  to  that 
of  a  hydrogen  bomb  (1010  kilowatt-hours)  (Battan  1961)  and  hurricanes  have 


rained  over  ninety-five  inches  over  a  spot  in  four  days  (Silver  Hill, 

Jamaica  in  November,  1909)  (Alaka  1968).  The  point  is  that  a  system  with 

this  much  energy  and  water  substance  is  no  local  accident,  but  could  be 

* 

appreciable  to  global  balances.  There  are  several  efforts  in  the  U.S. 
alone  on  global  atmospheric  simulation  on  high-speed,  large-storage 
digital  computers  for  both  short-term  weather  prediction  and  also  long¬ 
term  climatological  reconstruction  and  anticipated  evolution.1"  The  effort 
at  the  NOAA  Geophysical  Fluid  Dynamics  Laboratory  in  Princeton,  New  Jersey 
includes  a  program  under  Kurihara  (1971)  devoted  to  numerical  integration 
of  the  basic  equations  believed  to  model  a  tropical  cyclone.  Similarly, 
a  major  motivation  behind  the  TRW  tropical  cyclone  project  being  described 
here  is  its  possible  contribution  to  the  ARPA-sponsored  global  climatological 
analysis  under  way  at  Rand  Corporation  in  Santa  Monica,  California.  While 
there  are  other  familiar  rapidly  swirling  vortices  occurring  in  the  lower 
atmosphere  (tornadoes,  waterspouts,  dustdevils,  fi rewhirls,  and  swirls 
attending  volcanic  eruptions),  none  of  these  persist  long  enough  in  time  and 
occur  over  a  large  enough  area  to  warrant  consideration  for  incorporation 
in  current  feasible  global-scale  computations.  On  the  fastest  conventional 
computers  a  grid  no  finer  than  4°  x  5°  is  practical  for  the  Rand  program, 
and  even  on  an  advanced  computer  like  the  Illiac  IV  no  grid  finer  than 
1°  x  1°  seems  practical  (Rapp  1970). 

★ 

Because  hurricanes  have  radial  extents  of  several  hundred  miles  and  extend 
to  the  tropopause,  and  because  hurricanes  persist  for  weeks,  it  seems 
implausible  that  they  are  mere  accidents.  The  fact  that  they  occur  annually, 
mainly  in  the  autumn  after  the  long  summer  heating  of  the  tropical  oceans 
by  solar  radiation,  suggests  the  speculation  that  they  are  a  mechanism  for 
relaxing  energy  poleward  when  the  usual  Hadley-cell  mechanism  is  not 
sufficient.  Hurricanes  in  general  do  turn  poleward  after  drifting  westward 
in  the  trades. 

f0ne  effort  is  carried  on  at  the  Geophysical  Fluid  Dynamics  Laboratory  of 
NOAA  on  the  Forrestal  Campus  of  Princeton  University  in  New  Jersey;  this 
effort,  under  Dr.  J.  Smagorinsky,  has  Department  of  Commerce  funding 
(Smagorinsky  1963).  A  second  effort  is  conducted  by  Dr.  A.  Kasahara  under 
NS F  funding  at  Boulder,  Colorado  (Kasahara  and  Washington  1967).  The  third 
effort  is  supervised  by  Dr.  R.  Rapp  of  Rand  Corp.  in  Santa  Monica,  California 
under  ARPA  (Department  of  Defense)  funding;  this  effort  concerns  solution  of 
the  Mi ntz-Arakawa  model  on  the  Illiac  IV  computer  (Rapp  1970).  The  goals 
and  techniques  of  these  three  efforts  differ  appreciably. 
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I.B  SOME  OBSERVATIONAL  FACTS  ON  HURRICANES 

Every  year,  especially  in  the  late  summer  and  early  fall,  several  of 
the  many  very  large  depressions  in  the  trades  (specifically,  disturbances 
between  5r  and  15°  of  the  equator)  intensify  into  tropical  cyclones*  [winds 
in  excess  of  74  mph  by  convention]  over  warm  tropical  oceans  (usually  at 
least  26°C  to  27°C,  often  28°C  and  higher).  These  tropical  cyclones  are 
typically  a  thousand  miles  in  diameter  and  no  more  than  ten  miles  in  height ; 
they  often  persist  from  one  to  several  weeks,  traveling  westward  in  the 
barotropic  trades  at  about  twenty  miles  per  hours  before  moving  poleward 
at  greater  translational  speeds.  These  vortical  storms  are  cyclonic  in  the 
Northern  Hemisphere  and  anticyclonic  (North  Pole  reference)  in  the  Southern 
Hemisphere,  and  take  many  days  to  intensi fy--indicating  that  the  small 
rate  of  rotation  of  the  earth  is  the  source  of  angular  momentum  (Palmgn 
and  Newton  1969).  In  fact,  conservation  of  angular  momentum  in  itself 
indicates  that  a  fluid  particle  in  the  tropics  drawn  in  about  five  hundred 
miles  will  swirl  at  several  hundred  miles  per  hour.  Tropical  cyclones 
have  local  designations  around  the  globe  (e.g.,  hurricanes  in  the  North 
Atlantic,  typhoons  in  the  northwestern  Pacific,  papagallos  on  the  west  coast 
of  Central  America,  baguios  in  the  Philippine  Islands,  willy-willies  in 
Hustral ia ,  cyclones  in  the  northern  Indian  Ocean,  and  trovados  near 
Madagascar). 

Tropical  cyclones,  which  total  about  fifty  in  a  typical  year,  are 
known  over  all  oceans  except  the  South  Atlantic.  With  satellite  photo¬ 
graphy  the  inspection  of  broad  ocean  expanses  has  improved;  it  seems  that 
there  are  about  eight  tropical  storms  annually  in  the  North  Atlantic,  and 


"A  tropical  cydone  starts  out  as  a  tropical  disturbance  in  whirh  th*™ 

£1S£u?ation  and  P^haps  one  closed  isobar.  When  the  wind 
creases  to  about  2U  knots  and  there  is  more  than  one  closed  isobar 

to°more  than^fknots1  ^nd^h^  3  tropica1  feession.  When  the  wind  rises 
co  more  than  34  knots,  and  there  are  severa)  closed  isobars  it  becomps 

iMs  classTTT&rl?  TiT"1’  •  If  the  winds  excee<)  64  knots  <74  “i'es/hour) . 
location)"  (dITiSm!  Sf  ^  °r  (depend1n9 
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about  half  of  these  intensify  into  tropical  cyclones  (Meyer  1971).  Many 
Atlantic  hurricanes  can  be  traced  to  biweekly  disturbances  that  begin  as 
sandstorms  over  the  Sahara  (Palmen  and  Newton  1969);  these  cyclonic 
depressions  may  extend  1000  nautical  miles  and  drift  westward  at  up  to  ten 
knots.  Analysis  of  dust  samples  taken  on  Caribbean  isles  after  hurricane 
passage  reportedly  confirms  this  (Jennings  1970).  As  the  peak  of  the 
North  Atlantic  hurricane  season  approaches,  the  region  where  tropical  storms 
reach  hurricane  intensity  moves  eastward  from  the  Gul i‘  of  Mexico  and  the 
Caribbean  to  the  Cape  Verde  Islands;  as  the  hurricane  season  passes,  the 
spawning  ground  moves  westward  again  to  the  Caribbean  (Meyer  1971). 

It  is  known  that  hurricanes  form  where  there  is  sustained  local  convective 
activity  over  warm  tropical  seas  (so  that  air  lifted  on  a  moist  adiabat  remains 
warmer  than  the  undisturbed  ambient  up  to  12  km),  where  there  is  enhanced 
cyclonic  shear  (as  occurs  when  the  Intertropical  Convergence  Zone  lies  at  a 
considerable  distance  from  the  equator),  and  where  there  is  weak  vertical 
shear.  The  last  requirement  supposedly  explains  the  anomalous  cyclone 
season  for  the  northern  Indian  Ocean,  with  twin-peak  frequencies  of  occurrence 
(in  spring  and  fall)  with  a  relatively  uneventful  summer  season  (Palmtfn  and 
Newton  1969).  On  the  other  hand,  hurricanes  tend  to  weaken  over  land 
moderately  rapidly  [the  central  pressure  of  Camille  rose  from  905  to  990  mb 
in  about  thirteen  and  one-half  hours  after  land  fall  (Bradbury  1971)]. 

Tropical  cyclones  have  a  structure  characterized  in  the  mature  stage 
by  a  relatively  cloud-free  calm  (winds  usually  well  below  fifteen  mph)  eye 
of  about  ten  to  twenty  mile  radius.  The  eye  is  characterized  by  low 
pressure  at  sea  level  (often  below  960  mb)  and  high  temperatures  aloft 
(10°C  above  ambient).  The  eye  is  surrounded  by  an  eyewall.a  ten-mile  wide 


4. 

During  the  warmer  months,  at  least  one  easterly  wave  is  present  almost 
every  day  over  the  Atlantic.  In  that  region  an  average  of  only  eight 
disturbances  per  year  reach  tropical  storm  intensity  ...  and  about 
60  percent  of  these  achieve  hurricane  force  .  .  .  Thus  a  weak  disturbance 
Ml  a  poor  chance  of  becoming  a  tropical  storm,  but  one  tjjaQii  arhlowoJ 
trojjical  storm  intensity  has  an  excellent  prospect  of  becoming  a  full 
humcane1'  (Palmen  and  ftewton  1969,  p.  b'03).  There  appears  to  be  no 
criterion  such  that  once  a  developing  depression  exceeds  it,  it  will 
definitely  become  a  hurricane. 
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annulus  of  intense  convection,  torrential  rainfall,  and  deep,  thick 
cloudiness.  Outside  the  eye  are  convective  rainbands  that  appear  like 
p inwheels  or  logarithmic  spirals  in  some  satellite  photographs  and/or  radar 
displays  taken  from  above  the  storm.  The  principal  velocity  component  In 
much  of  the  storm  Is  azimuthal  (or  tangential);  the  vertical  velocity 
component  is  appreciable  In  the  eyewall.  There  Is  low-level  cyclonic 
inflow  and  (In  the  outer  regions)  high-level  antlcyclonlc  outflow  (relative 
to  an  observor  rotating  with  the  earth)  for  a  Northern  Hemisphere  tropical 
cyclone  (Palmen  and  Newton  1969).  While  there  Is  much  spray,  the  lowest 
few  hundred  feet  (at  least)  of  the  Inflow  layer  remain  cloud-free  In  as 
far  as  the  eyewall  (Rlehl  1954). 

It  Is  often  agreed  that  there  Is  slow  downward  motion  In  the  eye  and  In 
the  outer  regions  of  the  storm;  that  the  latent  heat  of  condensation 
reduces  the  density  In  the  eyewall  to  establish  a  large  radial  pressure 
deficit  relative  to  ambient  conditions,  from  hydrostatic  considerations; 
that  a  still  further  pressure  deficit  from  ambient  occurs  In  the  «ye  owing 
to  roughly  dry  adiabatic  recompression  of  air  that  has  risen  along  a  moist 
adlabat  In  the  eyewall;  and  that  a  cyclostrophlc  balance  (balance  of  radial 
pressure  gradient  and  centrifugal  force)  yields  a  good  estimate  of  the 
swirl  speeds. 

Cloud  cells  In  the  eyewall  are  typically  5  to  20  km  thick  and  the 
smaller  ones  can  rain  over  6  In./hr.  Spiral  bands  out  to  150  km  yield  0.4 
In/hr,  and  further  out,  0.1  In./hr-though  1  km  convective  cells  can  give 
much  heavier  rainfall.  In  the  low  rainfall  area  the  precipitation  Is 
probably  snow  that  turns  to  rain  at  the  melting  level  (Heyer  1971).  In 
addition  to  the  rainfall,  storm  surge,  and  large  waves  already  mentioned, 
tropical  cyclones  can  spawn  tornadoes  and  waterspouts  (Orton  1970).  Only 
near  the  center  are  tropical  cyclones  axlsymmetrlc  to  good  approximation; 
n«ar  the  outer  edges  there  is  asymmetry.  The  location  of  mo'*  '--»wned 
tornadoes,  maximum  rainfall,  the  most  Important  rainbands,  m..*  Highest 
winds  suggest  that,  for  North  Atlantic  hurricanes,  the  storm  Is  most  severe 
In  the  right  forward  quadrant,  with  respect  to  an  observor  looking  along  the 
direction  of  translation  (Hawkins  1971).  Clearly  the  additive  translational 
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contribution  to  the  azimuthal  wind*  it  the  plausible  explanation  of  the 
V  asymmetry  in  wind  speeds  (Rlehl  1954,  p.  290). 

itost  of  these  remarks  concern  the  nature  hurricane  and  its  dec  ay.  Hich 
of  what  is  written  about  intensification  is  tentative,  nebulous,  and 
labyrinthine.  One  point  about  which  there  is  obvious  indecision  is  the 
role  of  the  Intertropical  Convergence  Zone  (or  Trade  Confluence)  In  tropical 
cyclone  genesis.  Palmdn  and  Newton  (1969,  p.  503)  cl  to  one  source  (tan 
and  Miller  1969)  which  states  that  only  about  one-sixth  of  ".  .  Atlantic 
hurricanes  originate  as  perturbations  novlng  away  fro m  the  trade-confluence 
lone  in  the  Panama  region  .  .  while  another  source  (Palndn  and  Newton 
1969,  p.  476)  states  that  "the  great  najorlty  of  the  cyclones  of  the 
tropical  north  Pacific  fom  In  latitudes  south  of  6*N  but  intensify  rapidly 
only  in  the  tone  6*N  to  1S*N."  The  Own  and  Miller  boob,  in  de  enphasiiing 
the  role  of  the  doldruns,  is  at  variance  with  both  earlier  and  later  worti 
for  example,  Byers  (1944,  pp.  425-428)  notes  that  tropical  cyclogenesis 
occurs  only  in  those  locales  and  during  those  seasons  in  which  the  Inter¬ 
tropical  Convergence  Zone  (ITCZ)  has  been  displaced  far  from  the  equator— 
up  to  15*.  Often,  when  the  ITCZ  hos  been  displaced  this  far  froa  the 
equator,  below  about  6000  ft,  between  the  ITCZ  and  the  equator,  there  is 
a  westerly  flow  within  the  general  tropical  easterly  flow— a  possible 
source  of  enhanced  shear.  Byers  attributes  the  absence  of  hurricanes  In  the 
South  Atlantic  largely  to  the  failure  of  the  ITCZ  to  becoae  displaced  south 
of  the  equator,  even  in  February.  More  recently,  Charoey  (1971)  notes  that 
the  sane  conditions  which  are  important  to  tropical  cyclogenesis  (low-level 
convergence,  large  moisture  content  and  appreciable  Coriolis  force) 
and  the  sane  processes  active  In  topical  cyclogenesis  (conditional 
instability  of  the  second  Mnd--U$C)  are  also  Important  In  the  generation 
of  the  ITCZ.  The  CIS*  process  is  the  feeding  of  convective  activity  In  a 
swirling  flow  by  frirtionally  induced  Inflow  in  a  surface  boundary  layer; 
the  convection  so  sustained  by  the  low-level  moist  inflow  results  In  a 
radial  pressure  gradient  (through  local  lightening  of  the  air  by  condensa¬ 
tional  heat  'clease),  such  that  swirling  is  sustained  to  create  more 
inflow.  This  important  subject  will  be  addressed  again  later.  One  other 
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ffftiri  about  Intfns fflcatlon  Is  that  most  theories  propose  Ihii  *t>  tye  f* 

gradually  a  the  tropical  disturbance  grows  into  a  tropical  tyclone 
(Palmen  and  'teuton  1969). 

Pith  prediction  for  tropical -cyclones  Is  still  In  an  imperfect  state, 
prevailing  method  is  largely  historical-wHat  did  previous  similar  hurri- 
canes  *>a  in  similar  circumstances?  Carly  analytical  wort  sought  a  "’steering 
Invel*  in  the  ambient  winds  (Byers  (1944,  p.  447)  states  that  «...  the 
circulation  in  the  upper  el*  such  as  at  10,000  ft.,  determines  the  hurricane 
path  with  considerable  accuracy.”).  Later  a  -steering  layer”  concept  was 
found  more  satisfactory  (ftiehl  (1954,  p.  345)  states  that  ”...  tropical 
i\orm  npvfjn  |hp  iimijon  and  wljh  the  speed  of  the  s  teen  no  current. 

Uhlch  is  defined  as  the  pressured  weighted  mean  flew  from  the  surface  to 
***  nfc  over  a  band  8*  latitude  In  width  and  centered  on  the  storm”).  Still 
more  recent  *orfc  averages  over  the  depth  of  the  troposphere  from  100  to 
1000  mb  (Sanders  and  Burpee  1968).  There  era  many  special  circumstances; 
for  e»a*ple,  coexistent  binary  tropical  cyclones  In  the  s tm  hemisphere 
route  about  one  another  (rujluhara  effect)  (Brand  1970),  while  binary 
systems  In  different  hemispheres  tend  to  move  parallel  (Co*  and  Jag er 
1969)— as  suggested  by  classical  potential  theory  for  line  vortices. 

Prior  passage  of  a  previous  cyclone  can  also  have  an  effect  (Brand  1971). 
Tropical  cyclones  often  recurve  eastward  at  midlatitudes  along  the  western 
side  of  high-pressure  cells,  and  can  Interact  with  e*tr#tropical  cyclones 
(Palmdn  and  Newton  1969). 
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n.  NOOILS  Of  A  TROPICAL  CYCLONE 
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11. A  INTRODUCTION 

In  sorting  out  the  thermohydrodynamfcs  of  a  tropical  cyclone,  one  1» 
Uct4  with  understanding  the  Interaction  of  two  scales  of  phenomena.  One 
scale  Is  larger  cyclone  scale;  a  nature  Intense  tropical  cyclone  may 
easily  reach  radial  proportions  of  five-hundred  to  one-thousand  miles  before 
the  winds  subside  to  ambient.  The  smaller  cimulus  scale  Is  at  least  two, 
often  three,  orders  of  magnitude  smaller.  The  cyclone  must  feed  the 
cumulus  scale,  which  In  turn  sustains  the  cyclone  scale.  In  a  cooperative 
Interdependence. 

Before  turning  to  this  Interaction  problem,  It  Is  probably  worthwhile 
to  categorlte  most  current  research  In  terms  of  these  two  scales.  Almost 
all  current  theoretical  research  on  tropical  cyclones  concentrates  on 
perameterltlng  the  cumilus  convection;  only  by  so  doing  can  the  hurricane 
be  properly  described  on  a  high-speed  digital  computer.  According  to  most 
workers,  the  gross  cyclone-scale  themohydrodynvHcs  has  already  been 
essentially  and  correctly  outlined  by  Rlehl  and  Malkus  In  a  series  of 
articles  and  books  In  the  late  1950's  and  early  1960's  (Rlehl  1954;  Malkus 
195b;  Malkus  and  Rlehl  1960).  A  divergent  point  of  view  about  cyclone- 
scale  thcmohydrodyneml cs  has  been  set  forth  by  Carrier  and  his  co-workers 
In  a  series  of  articles  published  within  the  last  two  years  (Carrier  1970; 
Dergarafcedian  and  fer.dell  1970;  Dergarabedlan  and  Fendell  1971;  Carrier. 
Hammond,  and  Georg®  1971;  Carrier  1971a;  Carrier  1971b).  These  articles 
delineate  overall  dynamics  and  thermodynamics,  and  Imply  (In  contrast  to 
Rfehl-ttalkus)  no  major  augmentation  of  ambient  heat/moisture  transfer 
from  the  ocean  is  needed  to  eaplaln  hurricanes. 

The  controversy  between  the  Rlehl-Nalkus  and  Carrier  theories  about 
cyclone-scale  thereohydrodynamics  has  been  briefly  discussed  before,  but 
shall  be  reviewed  here  for  two  reasons,  first,  recent  publications  have 
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clarified  the  disagreements  particularly  clearly.  Second  --  and  most 
Important  —  there  seems  a  strong  possibility  that  the  same  physical 
processes  operating  on  the  cyclone  scale  are  also  operating  (properly 
scaled  down)  on  the  cumulus  scale.  This  makes  it  all  the  more  important 
to  appreciate  the  structure  and  processes  operative  in  a  tropical  cyclone 
in  the  large. 

II. B  THE  CARRIER  MODEL 

In  this  section, mainly  the  structure  of  a  mature,  fully  developed 
tropical  cyclone,  taken  to  be  adequately  modeled  as  quasi  steady  and  axi- 
synmetric,  will  be  studied.  In  later  sections  transient  analysis  aimed 
at  establishing  how  such  a  severe  vortical  storm  is  established  will  be 
undertaken. 

The  Carrier  model,  on  the  basis  of  subdividing  the  tropical  cyclone 
into  segments  where  different  processes  and  scales  predominate,  is  a  four- 
part  analysis.  The  four  regions,  Indicated  in  Figure  II.B.l,  are  the 
throughput  supply  I,  the  frictional  boundary  layer  II,  the  eyewall  and 
efflux  region  III,  and  the  eye  IV.  Some  of  this  subdivision  is  conventional, 
some  not.  Besides  clarifying  locally  dominant  physical  processes.  It 
permits  retention  of  the  minimal  number  of  terms  in  locally  valid  quanti¬ 
tative  formulation;  this  procedure  simplifies  the  mathematical  solution  in 
a  manner  unavailable  to  any  direct  finite-differencing  of  uniformly  valid 
equations. 

The  Carrier  model  is  closed  for  convenience  —  there  is  no  very 
significant  amount  of  mass  convected  across  any  boundary.  The  cylindrical- 
like  volume  encompassing  the  entire  storm  has  the  sea  surface  for  its 
bottom;  Its  sides  lie  far  enough  from  the  center  (about  500  to  1000  mi)  so 
that  the  winds  are  virtually  reduced  to  ambient,  and  the  swirl  relative 
to  the  earth  is  taken  as  zero.  The  top  of  the  storm  is  taken  to  be  that 
height  at  which  sea-level  air  In  the  outer  part  of  the  storm,  if  lifted 
rapidly  so  that  the  total  enthalpy  of  a  fluid  particle  remained  constant 
because  relatively  slow  ambient-maintaining  processes  would  not  have  time 
to  act,  would  no  longer  be  unstable  relative  to  the  local  ambient  air.  Such 
sea-level  air  would,  of  course,  rise  dry  adiabatically  until  saturated; 
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Fig.  II.B.l  This  conjectured  configuration  of  a  mature  hurricane  with 
rough  order-of-magni tude  dimensions  is  not  drawn  to  scale.  The  subdomains 
are:  I,  throughput  supply,  a  region  of  rapid  swirl  and  slow  downdraft; 

II,  frictional  boundary  layer;  III,  eyewall;  and  IV,  eye.  Across  the 
boundary  layer  II  there  is  about  a  100  mb  drop  and  across  I,  a  further 
200  mb  drop;  at  the  top  of  the  hurricane  the  pressure  Is  about  150  mb,  l.e. , 
the  top  is  near  the  tropopause. 
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thereafter,  enough  precipitation  would  fall  out  to  leave  the  air  Just 

; 

saturated  at  the  local  temperature  and  pressure,  the  air  retaining  the 
latent  heat  of  phase  transition.  Such  a  Incus  of  states  Is  conventionally 
referred  to  as  the  moist  adiabat.  The  vertical  profile  of  the  equivalent 
potential  temperature  in  a  typical  autumnal  tropical  ambient  Is  such  that 
it  decreases  with  increasing  height  to  roughly  330°K  at  650  mb,  then 
increases  and  recovers  Its  sea-level  value  of  about  350#K  at  about  150  N>. 
The  autumnal  tropical  ambient  Is  conditionally  unstable  in  that  a  particle 
at  any  height  lifted  dry  adlabatlcally  attains  a  potential  temperature 
lower  than  that  of  the  tropical  ambient  and  hence  would  return  to  Its 
Initial  position;  but  any  sea-level  particle  displaced  vertically  enough 
for  the  onset  of  condensation  would  continue  to  rise  many  Mlweters.  This 
"Instability  lid"  lies  at  so  great  a  height  that  there  Is  virtually 
negligible  swirl,  as  explained  below;  the  ambient  pressure  and  temperature 
at  this  height  Is  taken  to  describe  all  radial  positions  at  this  height, 
from  the  center  to  the  outer  edge.  Thus,  the  top  of  the  storm  Is  an  Iso¬ 
thermal,  Isobarlc  lid  with  no  water  vapor  content  for  current  purposes. 

The  relatively  slow  processes  that  maintain  the  tropical  ambient,  to 
which  allusion  has  just  been  made,  are  cumulus  convection,  turbulent 
mixing,  and  radlatlonal  transfer.  While  these  are  quail Utlvely  easy  to 
describe,  precise  quantitative  formulation  Is  very  formidable. 

Discussion  now  turns  to  describing  each  of  the  four  regions  comprising 
the  tropical  cyclone  In  some  detail. 

In  region  I  there  Is  warm  moist  air  typical  In  stratification  of  the 
ambient  atmosphere  in  which  the  hurricane  was  generated.  This  air  spun  up 
under  conservation  of  angular  momentum  as  It  moved  In  toward  the  axis  of 
symmetry  during  the  formative  stage.  As  the  mature  stage  was  approached, 
a  gradient-wind  balance  of  pressure,  Coriolis,  and  centrifugal  forces 
choked  off  any  further  Inflow;  the  Inflow  Is  only  enough  to  prevent  the 
eyewall  III  from  diffusing  outward,  and  that  requires  cnly  an  exceedingly 
small  radial  flow.  The  air  In  I,  then.  Is  rapidly  swirling,  the  azimuthal 
velocity  component  greatly  Increasing  and  the  pressure  greatly  decreasing 
from  the  e'je  to  the  center.  Under  such  a  radial  profile  for  the  swirl. 
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lhff*  I*  a  smalt  dewnflu*  fro*  the  th'Ougnput  supply  I  ln|n  the  frictional 
bou«MJiry  layer  ||.  The  smell  downf lu*  leads  to  a  large  not  mm  flu*  Into 
thf  friction* I  layer  tftmi  of  the  iff*  involved,  furthermore,  the  down* 
draft  Is  only  4  9*011  average  because  locally  4*d  transiently  there  Is 
intense  elective  activity  by  which  clouds  form  end  rein  fells.  Ihe 
clouds  ere  strained  by  the  rapid  swirl  Into  the  spiral  bends  s*en  on  radar 
screens  or  In  satellite  photographs.  The  spiral  bends  9ive  vlsuellietfon 
to  parts  of  ihe  strain  pattern,  rather  th*r  streamline  pattern. 

In  the  frictional  boundary  layer  It,  the  only  region  In  which  angular 
nomentww  Is  not  conserved  but  Is  partially  lost  to  the  sea,  there  Is 
appreciable  Influ*.  In  fact,  the  atimuthel  and  radial  velocity  components 
are  of  comparable  magnitude;  typically,  for  flaed  radial  poslti&n,  the 
rating  inflow  speed  a*,  apy  a*ial  position  in  the  boundary  layer  Is  about 
one-thlrd  the  mvInwR  at  muthel  speed,  [this  fraction  Is  about  the  one 
reported  by  Hughes  (I9S?)  from  flight  penetration  of  hurricanes  at  altitudes 
of  iOOn  feet  or  less.]  The  vertical  velocity  component  Is  ouch  smaller. 

The  reason  for  the  Inflow  Is,  of  course,  that  the  no-slip  boundary  condition 
reduces  the  centrifugal  acceleration,  and  a  relatively  uncompensated  pressure 
gradient  drives  the  fluid  tmard  the  a*ls  of  sy*netry  (so-called  *tee  cup 
effect*).  Ter  from  the  e«1i  In  II  the  classical  balance  of  the  linear 
Clman  layer  (friction,  pressure,  and  Coriolis  fortes)  suffices;  since  the 
d&ndreft  from  I  to  II  is  probably  fairly  Independent  of  radial  position 
(especially  far  from  the  eyewall)  for  twirl  distributions  of  practical 
Interest,  three-quarters  of  the  flu*  Inward  through  It  cones  from  downflu* 
across  the  Interface  between  It  and  1  where  (rJ2)  <  r  e  r0,  where  rQ  Is 
the  radial  e*lcnt  of  the  storm.  Closer  In  to  the  avis  the  nonlinear 
accelerations,  especially  the  radial  acceleration  of  radial  and  tangential 
remen t’^  and  centrifugal  acceleration,  must  enter. 

M  the  pressure  gradient  in  I  lets  more  air  slnl  Into  II,  the  Influ* 
in  II  drives  the  boundary  layer  air  moderately  rapidly  up  a  cloudy  eyewall 
HI.  In  ?*e  eyewall  hydrostatic  and  eye  I  os  trophic  approbations  hold; 

!*e  locus  0*  tN?f*©dyne*ic  state?  Is  the  moist  adfabat  based  on  sea-level 
conditions  in  the  eyewall.  The  swirl  near  the  top  of  til  is  so  reduced 


thn,  in  tht  outflow,  the  air  seems  to  an  observor  on  earth  to  be  rotating 
opposite  In  sense  to  the  rotation  In  I  and  II  (which  Is  cyclonic  In  the 
Northern  Hemisphere  and  antlcyclonfc  In  the  Southern).  The  air  In  III 
slips  over  the  air  In  I  with  no  Interaction;  there  Is  no  large  radial 
pressure  gradient  In  much  of  111,  unlike  I. 

At  low  altitudes  the  eyewall  flushes  moist  air  out  of  IV,  and  at 
high  altitude  ralned-out  air  is  entrained  into  IV.  In  time  the  eye  becomes 
better  defined;  It  Is  the  central  core  In  which  relatively  dry  air  sinks, 

Is  warned  by  compression,  and  is  entrained  out  Into  the  eyewall  or  re¬ 
circulated  within  the  eye  (stagnation  at  the  base  of  the  eye  would  permit 
transport  processes  to  cool  the  eye).  The  relatively  light  eye  permits 
much  greater  pressure  deficits  from  ambient  and  hence  supports  much  higher 
swirling  speeds.  At  a  fixed  altitude,  the  density  In  the  eye  is  less  than 
that  In  the  eyewall,  which  In  turn  Is  less  than  that  of  the  m*1ent  gas  at 
the  storm  edge.  Since  the  hydrostatic  approximation  Is  uniformly  valid, 
spatially  and  temporally.  In  a  hurricane,  the  sea-surface  eye  pressure  is 
less  than  the  see-surface  eyewall  pressure,  which  is  less  than  the  sea- 
surface  ambient  pressure. 

Carrier's  model  thus  pictures  the  tropical  cyclone  as  a  once-th rough 
process  In  idilch  a  -fuel  supply"  —  the  warm  moist  air  In  I,  a  part  of  the 
tropical  cyclone  at  Its  Inception  and  convected  with  the  storm  -  Is  slowly 
exhausted.  The  storm  weakens  because  the  air  drifting  down  Into  the 
frictional  lqyer  tward  the  end  Is  typical  of  the  higher  tropical  environ¬ 
ment  and  hence  of  lower  equivalent  potential  Umperature.  Eventually  the 
fuel  supply  Is  exhausted,  and  the  boundary  between  III  and  I  sinks  toward 
II.  Some  models  (e.g.,  £1  lessen  and  Klelnschntdt  1557)  picture  a  recirc¬ 
ulation  through  the  storm  of  outflow  air;  the  storm  does  not  survive  long 
enoigh  for  this,  nor  could  such  air  maintain  the  storm. 

There  Is  one  Important  omission  to  the  foregoing  description  that  has 
been  Intentionally  deferred:  the  energetics  of  the  surface  frictional 
layer  and  associated  questions  of  air/sea  transfer  of  latent  and  sensible 
heat.  The  relevant  quantity  to  consider  is  the  total  stagnation 
enthalpy  (the  sums  pf  static  enthalpy,  the  heat  associated  with  condensible 
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moisture,  gravitational  potential  energy,  and  kinetic  energy  contributions). 
This  quantity  (a  generalization  of  the  equivalent  potential  temperature) 
is  conserved  at  roughly  its  ambient  stratification  throughout  regions  I  and 
II;  therefore,  it  Is  described  by  a  profile  that  decreases  with  height  from 
1000  mb  to  650  mb,  and  then  Increases  with  height,  as  mentioned  above.  The 
implication  is  that  the  heat  and  mass  transfer  from  the  ocean  to  the 
atmosphere  is  about  the  same  within  the  hurricane  as  In  the  ambient.  This 
transfer  helps  compensate  for  the  rain-out  in  the  spiral  bands  and  helps 
maintain  the  warm,  moist  nature  of  the  air  in  I.  [Occasionally  the  Carrier 
model  is  still  grossly  misrepresented  as  proposing  adiabatic  conditions 
(constant  total  stagnation  enthalpy  In  II  so  the  net  heat  and  mass  transfer 
from  sea  to  air  is  zero);  such  a  solution  cannot  possibly  satisfy  the 
parabolic  boundary-value  problem  describing  the  energetics  of  the  frictional 
boundary  layer  because  It  obviously  violates  the  Initial  condition  at 
r  =  rQ,  the  outer  edge.  In  fact,  if  the  supplemental  flux  from  the  ocean 
is  entirely  eliminated,  as  from  passage  over  land,  the  spin-down  time  is 
0(a/v  /2fi1/2)  where  the  edc(y  viscosity  v  =  10‘2  m12/hr,  the  normal  component 
of  the  rotation  of  the  earth  n  i  2  X  10"1  hr-1,  and  the  height  of  the 
throughput  supply  a  *  1  ml  —  so  the  spin-down  time  Is  half  a  day  to  two 
days.]  The  model  of  total  stagnation  enthalpy  fixed  at  its  ambient 
stratification  breaks  down  In  the  eyewall  III;  there  the  vertical  velocity 
component  is  at  least  one,  probably  two  orders  of  magnitude  larger  than 
the  relatively  small  downdrlft  Into  the  boundary  layer;  the  result  Is  that 
convection  dominates  the  slow  ambient-sustaining  processes  so  the  total 
stagnation  enthalpy  Is  virtually  constant  at  Its  sea-level  value,  which  Is 
roughly  its  ambient  sea-level  value. 

Three  specific  analyses  carried  out  by  Carrier  and  his  co-workers  to 
corroborate  aspects  of  this  quasi  steady  model  of  the  mature  tropical 
cyclone  are  now  briefly  reviewed.  These  Involve  maximum  swirl  speed 
estimation,  the  dynamics  of  a  nonlinear  Ekman  layer,  and  the  energetics  of 
the  surface  frictional  layer.  These  problems  are  not  difficult  to  formulate 
nor,  for  the  accuracy  of  result  required,  are  they  difficult  to  solve.  Much 
novel  and  valuable  information  about  tropical  cyclones  is  attainable  without 
large-scale  computation. 
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II.B.l  MAXIMUM  SWIRL  SPEED  ESTIMATE 

An  upper  and  lower  bound  on  the  central  pressure  deficit  achievable 
in  a  known  spawning  atmosphere  will  now  be  set  forth  by  use  of  the  hurricane 
model  just  presented,  of  hydrostatics,  and  of  the  thermodynamics  of  moist 
and  dry  air.  Specifically,  the  weights  of  various  columns  of  air  in  the 
storm  will  be  determined  in  light  of  different  moisture  content  and  thermo¬ 
dynamic  processes  involved.  The  bounds  on  the  central  pressure  deficit  can 
then  be  translated  into  an  estimate  of  bounds  on  the  maximum  swirl  speed 
through  dynamics  (the  radial  momentum  equation).  Fletcher  (1955)  had 
suggested  use  of  the  cyclostrophic  balance  once  pressure  deficits  were 
known,  and  Malkus  (1968)  had  suggested  that  pressure  deficits  could  be 
calculated  from  moi^t  adiabatic  considerations  for  the  eyewall  and  dry 
adiabatic  considerations  for  the  eye.  Here  the  concepts  are  combined  to 
achieve  quantitative  bounds,  but  just  as  important,  to  demonstrate  that 
hurricane  speeds  could  be  achieved  without  requiring  any  augmenting  enthalpy 
transfer  from  the  ocean  whatever. 

The  first  step  is  to  neglect  the  frictional  boundary  layer  II,  which 
is  relatively  thin  and  across  which,  except  for  hydrostatic  variations, 
the  pressure  does  not  change  according  to  lowest-order  boundary  layer  theory. 

The  variation  of  pressure  p,  density  p,  and  temperature  T  with  height 
above  the  ocean  z,  for  any  ambient  tropical  atmosphere  in  which  a  hurricane 
forms,  may  be  computed  from 


pa  =  pa  Ra  T 

(a  =  dry  air); 

(II.B.l) 

Pv  *  pv  Ra  T/o 

(v  =  water  vapor;  o  =  0.622); 

(II. B. 2) 

P  =  Pa  +  Pv,  p 

=  pa  +  V  pv  =  P(T)  (RH^ 

(II. B. 3) 

d£ 

3z  ■  p  s; 

(II. B. 4) 

T  *  f(p),  RH  = 

g(p). 

(II. B. 5) 
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where  the  temperature  profile  f(p)  and  the  relative  humidity  (RH)  profile 
g(p)  are  taken  as  known  from  measurement.  The  saturation  pressure  P(T)  is 
wel 1- tabulated  for  vapor  and  liquid  phases  above  freezing,  and  vapor  and 
solid  below  freezing  (Keenan  and  Keyes  1936);  a  convenient  and  accurate 
expression  for  P(T)  in  mb  was  given  by  Tetens  (Murray  1967): 


P(T)  =  6.1078  exp 


a  (T  -  273.16) 
(T  -  b) 


( 1 1 .  B .  6a ) 


a  =  21.8745584  )  a  =  17.2693882  ) 

>  over  ice;  >  over  water,  (II.B.6b) 

b  =  7.66  )  b  =  35.86 


where  T  is  ,  °K.  The  integration  proceeds  from  the  sea  level  upward  in 

altitude  z;  data  typically  extend  from  about  1000  mb  to  150  mb.  The  top 
of  the  storm  is  normally  taken  as  the  height  at  which  the  ambient  total 
stagnation  enthalpy  (for  which  the  kinetic  energy  contribution  is  negligible) 
recovers  its  sea-level  value,  as  noted  earlier;  here,  however,  a  slightly 
different  procedure  explained  below  will  be  used. 

In  a  fully  developed  storm  the  air  rising  up  the  eyewall  in  Carrier's 
model  follows  a  moist  adiabat  based  on  the  sea-level  ambient  state  (until  . 

late  in  the  storm  when  an  ambient  state  above  sea-level  should  serve  as  1 

the  reference  state  for  the  moist  adiabat,  but  by  then  the  storm  has 
weakened  from  the  maximum  intensity  level  of  interest  here).  Thus  for  the 
eyewall  one  integrates  dH  =  0  where  the  total  stagnation  enthalpy  H  is 
given  by 


H  =  cpT  +  LY  +  gz  +  q2/2  (II. B. 7) 

where  q  is  the  wind  speed  relative  to  the  earth,  L  the  latent  heat  of 
condensation  per  unit  mass,  and  Y  the  mass  fraction  of  water  vapor  (p  /p). 
After  manipulation  with  (I I .B.l )-( II  .B.4) ,  dH  =  0  may  be  written  (L  held 
constant) 


dT 

dp 


c 

"  px2dl 


(II. B. 8) 
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where  the  [d(qV2)/dp]  contribution  Is  negligible  and 


x  =  1  .  iLr.  .q)p  =  1  -  p  -  pa  .  (II.B..M 

A  slightly  different  mode  of  derivation  gives  a  similar  result: 


dT 

3p 


v 


Lo 
P  P, 


a  "  "a 
fo  '+  ft  T 


CP  + 


dP 

37 


( II cB.10) 


This  alternate  form  Is  used  in  any  calculations  reported  here.  The 
procedure  Is  to  use  the  dry  adiabatic  relation  T  ^  pY^Y  "  1  \  Y  .  i.4$ 
from  sea-level  ambient  conditions  [T(ps)  ■  T  ,  RH(p$)  -  RH$,  where 
subscript  s  denotes  sea-level  ambient  conditions  (given)];  where  RH  ■  1, 
one  switches  to  the  moist  adlabat  (II.B.10)  and  continues.  The  Integration 
Is  terminated  at  that  pressure  p  for  which  the  temperature  calculated  from 
the  moist  adlabat  (II.B.10)  and  from  the  ambient  (II.B.1)-(II.B.5)  are 
equal;  this  temperature  Is  denoted  Tj,  and  the  height  above  sea-level 
at  which  T1  occurs  Is  denoted  z}  (the  H11dM  on  the  cyclone).  One  then 
Integrates  (II.B.1)-(II.B.4)  and  (II.B.10)  from  z  -  z,  (where  T  -  T,  and 
P  =  Pj)  to  z  «  0;  RH  «  1  during  this  Integration  since  typically  sea-level 
autumnal  tropical  air  becomes  saturated  when  raised  even  20  mb.  If  no  eye 
existed  In  the  vortex  —  as  seems  to  be  the  case  for  some  tornadoes  and 
waterspouts  —  then  the  just-calculated  p(z  ■  0)  =  pg,  p(z  ■  0)  =  p  would 
characterize  conditions  at  the  center  of  the  vortex. 

In  a  mature  hurricane  a  pressure  deficit  In  excess  of  (p  -  p  )  Is 
achieved  by  having  rain-out  air  entrained  from  the  eyewall  sink  1n*a 
relatively  dry  eye  under  adiabatic  recompression.  Thus  In  a  hurricane 
(Ps  ■  Pe)  a  lower  bound  on  the  central  pressure  deficit.  For  an  upper 

bound  on  the  deficit  that  may  be  achieved,  one  may  adopt  the  Idealized  model 
that  the  eye  Is  completely  dry  (so  no  compresslonal  heat  Is  lost  to  re¬ 
evaporation)  and  that  the  air  entrained  Into  the  eye  Is  drawn  from  the  top 
of  the  eyewall  (or.  In  a ny  case,  has  T  •  Tr  p  -  pv  y  -  o  at  z  ■  z,). 
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The  relevant  equations  are  (II.B.1)-(II.B.4),  RH  «  0,  and  T  %  p**  *  D/t; 
integration  in  the  direction  of  decreasing  i  yields  p(z  >  0)  s  p  (<p  <  p  ) 
and  t (2  ■  0)  •  oc(<of  <  o%)  ~  the  density  discrepancies  so  calculated  are 
at  most  twenty- five  percent  and  Fletcher  (1955)  estimates  the  density  does 
not  vary  by  even  fifteen  percent,  so  the  density  nay  be  held  constant  at 
Its  ambient  value  throughout  the  tynanlcal  calculations  now  discussed. 

If  one  adopts  the  cyclostrophfc  balance,  holds  0  constant  at  (say) 

*V  and  (since  the  core  Is  observed  not  to  rotate)  lets 


v(r)  • 


0  0  <  r  <  R 

(v)wx  lP/r)"  R  <  r  <  • 


(II.B.11) 


then 


(II.B.12) 


First,  for  a  one-cell  vortex  (when  there  Is  no  aye  so  the  nolst  adlabat 
calculation  Is  appropriate  all  the  way  to  the  axis),  a  rlgld-boty-llke 
rotation  lies  near  the  core  so  then 


v(r) 


*v*max  fr/R) 

0  <  r  <  R 

*v,max  *R/r,n 

R  <  r  <  • 

(II. 6. 13) 


and  from  the  cyclostrophlc  balance 


(II.B.14) 


Next,  although  power-law  decays  of  swirl  with  radial  distance  are  frequently 
adopted  and  suffice  for  current  purposes,  it  will  become  evident  that  other 
forms  are  at  least  as  plausible,  and  more  convenient,  for  r  >  R.  In  any 
case.  Miliar  (1967)  suggested  from  limited  data  that  0.5  «  n~<  0.65  usually 
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suffices;  earlier,  Byers  (1944,  p.  435)  recommended  n  •  0.5  and  Hughes 
(1952)  considered  n  »  0.6  -0.7  suitable  for  an  average  hurricane.  The 

upshot  is  that  estimates  will  be  wade  on  the  limits  n  *  0.5  and  n  *  1.0. 

rinally,  the  gradient  wind  equation  would  probably  be  more  appropriate 
than  the  cyclostroph'c  equation,  but  the  more  complicated  formula  would 
give  maximum  speeds  reduced  by  only  five  percent  from  those  obtained  fro* 
the  simple  forms  (II.B.12)  and  (It .B. 14) . 

Tor  the  typical  ambient  West  Indies  tephigram  for  September  given  by 
Jordan  (1957)  (extending  to  the  130  mb  level),  one  calculates  a  pressure 
deficit  of  58  mb  for  a  one-cell  storm  (no  *ye)  and  137.5  mb  for  a  two-cell 
storm  (with  an  eye).  These  deficits  translate  Into  bounds  on  (v)^  of 
130.2  mph  for  n  *  0.5  and  159.5  mph  for  n  ■  1.0  for  the  storm  without  an  eye, 

and  224.8  mph  for  n  «  0.5  and  346  mph  for  n  *  1.0  for  a  storm  with  an  eye. 

further  results  are  given  In  Figures  (II.B.2)-(I1 .B.4). 

I  *  .82  THE  SMI  Rl-DI  VERCENCE  RELATION  FOR  THE  FRICTIONAL  BOUNDARY  LAYER 

Because  the  maxlmm  speed  achieved  in  a  tropical  cyclone  Is  rarely 
much  over  200  mph,  as  demonstrated  in  Section  II.B.l,  the  Hach  mrniber  rarely 
reaches  even  0.3.  Hence,  when  examining  the  dynamics  (as  opposed  to  the 
energetics),  an  Incompressible  constant-property  model  suffices. 

A  steady  axlsynmetrlc  fl<w  of  an  incompressible  fluid  is  now  studied 
to  confirm  the  crucial  point  that,  under  rapid  swirling,  there  is  downflux 
from  region  1  to  region  11,  and  sufficient  downflux  enters  the  surface 
frictional  layer  to  account  for  the  mass  flux  up  the  eyewall.  The  analysis 
will  be  carried  out  in  a  noninertial  coordinate  system  rotating  at  the 
constant  speed  of  that  component  of  the  rotation  of  thu  earth  which  is 
normal  to  the  local  tangent  plane.  Because  the  boundary-layer  divergence 
under  an  impressed  swirl  (the  major  constraint  furnished  by  the  boundary 
layer  on  the  inviscid  flow  above  it)  is  relatively  small  In  magnitude, 
careful  formulation  and  solution  of  the  coupled  quasll Inear  parabolic 
partial  differential  equations  and  boundary  conditions  describing  the  layer 
is  required. 
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Fig.  11.B.2  This  tephlgram  gives  thermodynamic  loci  determining 
upper  and  lower  bounds  on  sea-level  central  pressure  deficits  In  a  hurricane 
spawned  In  a  known  ambient  environment.  The  ambient  pressure-temperature 
curve  Is  based  on  data  For  the  Caribbean  In  September  by  Jordan  (1957). 

The  curve  labeled  moist  adlabat  Is  based  on  having  sea-level  a*61ent  air 
rise  dry  adiabatlcal ly  until  saturation,  and  thenceforth  moist  adlabatlcally. 
The  sea-level  pressure  associated  with  such  a  column  of  gas  gives  a  lower 
bound  on  the  central  pressure  deficit  from  ambient.  An  upper  bound  on  the 
deficit  is  furnished  by  having  the  air  that  rose  on  the  so-called  moist 
adlabat,  recompressed  dry  adiabatlcal Ty  from  the  top  of  the  storm  back 
down  to  sea  level.  Altitudes  are  associated  with  the  thermodynamic  state 
by  use  of  Hydrostatics  and  the  equations  of  state  for  dry  air  and  water 
vapor. 
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Fig.  II. B. 3  The  pressure-altitude  relation  for  each  of  the  three 
thermodynamic  loci  of  Fig.  II. B. 2  is  given  here  more  explicitly.  Clearly 
a  two-cell  storm  with  a  dry  adiabatic  core  can  support  much  higher  swirling 
speeds  than  a  one-cell  storm,  in  which  the  so-called  moist  adiabatic 
ascent  characterizes  the  flow  at  the  axis  of  rotation. 


The  relevant  equations  are 


v  *  q  B  o  ,  (11.0.15) 

v(q2/2)  ►  (7  x  q)  x  q  +  2?^  x  q  -  -vjl  -  v  Vx(?  x  q) ,  (II. 0.16) 

where  p  =  (p/p)  +  (ne  x  r)  /2  ♦  gz,  the  gravitational  acceleration  g  ■  -gz, 
the  velocity  In  ncninertlal  coordinates  v  =  n  x  r  ♦  q,  the  component  of 
the  rotation  of  the  earth  normal  to  the  local  tangent  plane  n  ■  nzt  and 
the  kinematic  viscosity  (later  given  eddy-diffusivity  values)  is  v. 

Nondimensionalization  Is  effected  by  letting  q'  •  q/(t  n)^2, 

P1  =  P/(foq)»  r'  =  r/('YQ/n)^^%  and  E  *  v/fQ  where  the  Ekman  number  E  <<  1 
and  yo  characterizes  the  circulation  away  from  the  boundary  (such  as  the 
maximum  swirl  speed  times  the  radius  at  which  It  occurs).  Dropping  primes, 
one  has 


*  •  Q  •  (II. 0.17) 

v(q  /2)  +  (vx  q)  x  q  ♦  2z  x  q  =  -vp  -  E  Vx(v  x  q).  (II. 0.18) 

These  equations  are  studied  In  axlsymmetrlc  cylindrical  polar  coordinates 

q  =  ur+ve+wz,  r-rr  +  zz  .  (II. 8. 19) 

Away  from  the  boundary  (l.e.,  In  region  I)  the  following  expansions 
are  ad?oted: 


P  =  i»(r,z)  +  ....  v  =  V(r,z)  ♦  ....  w  =  E1/2  W(r,z)  +  ...»  u  «  o(E1/2). 

(II .0.20) 

Substitution  of  (II . B. 20)  in  (II. 0.17)  and  (I .0.18)  gives  the  gradient- 
wind  equation: 


B  o,  W2  »  0,  *r  -  2V  ♦  V2/r  . 


(II. 8.21) 
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Subscripts  r  and  2  her#  («nd  x  «nd  c  below)  denote  partial  differentiation. 
The  axiaMy  Invariant  solution  Is  complete  when  Hr)  or  V(r)  Is  specified 
[here  V(r)  will  be  given);  M(r)  Is  found  by  etching  the  solution  to 
(M.B.2I)  to  the  solution  for  the  frictional  leyer  II,  end  In  this  sense 
W(r)  is  determined  by  the  boundary-layer  dynamics. 

If  c  •  2  r  [which  Implies  that  the  frictional  Uyer  Is  0(E1/2)  in 
thickness]  end  If  near  the  boundary 


U  •  ub(r.O  •  ....  ,  .  „b(r.<)  ♦  ....  „  .  c'/V(r.()  *  ....  p  .  p.(r.c) . 


(11.8.22) 


then  the  axial  component  of  the  momentum  conservation  equation  degenerates 
to  (*Pb*H)  *  0  in  conventional  fashion,  so  the  pressure  field  in  the 
boundary  layer  is  known  from  (II. 8. 21).  If 


<  *  rvb,  v  •  r  V.  «  •  rub,  x  •  r2,  w  •  2’]/^b.  £  •  2,/2  C,  (II. 8. 23) 

then  in  terns  of  dimensional  quantities  i  •  ru/r0,  e  *  rv/v  ,  w  •  w/(2ftv)^2, 
t  •  2/(v/2(r)1/2,  X  *  Qr2/f0.  the  boundary  layer  thickness  is  0(v/o),/2. 

In  terms  of  quantities  introduced  In  (II. 8.23),  one  has  frwi  (II. 8.17)  aid 
(II. 8. 18) 


<x*wC*°*  (II. 8. 24) 

::x  ♦  wtc  ♦  (v2  -  *2  -  42)/2x  -  (,  -  v)  -  .  0,  (II. 8.25) 

*♦*  ■  0  •  (II. 8. 26) 

hatching  of  expansions  gives 

:-0,  *  •  r(x)  given,  (II. 8. 27) 

and  at  *  «  0  no-slin  conditions  are  adopted: 

4  '  0:  *  *  "  *  *  *  0  .  (II. 8. 28) 
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Initial  conditions  arc  conveniently  given  by  noting  that  at  x  ■  x  ,  for  x 
large  enough,  solution  Is  given  by  discarding  all  the  nonlinear  terns  and° 
retaining  the  linear  eolations  treated  by  Elman,  In  which  x  enters  para¬ 
metrically  only.  The  solution  to  the  balance  of  Coriolis,  pressure,  and 
friction  forces  Is  well  known: 

*  •  -  [t(x)3  sin  (2*,/2c)  exp  (-2‘1/2c),  (II.B.29) 

♦  ■  (*(*)]  [1  -  cos  (2‘,/20  exp  (-2*,/2c)]  ,  (II .B . 30) 

w  *  2'1/2  (vx(x)]  {l  -  [s1n(2‘1/20  *  cos(2‘1/20]  [exp(-2‘1/20]  }  . 

(II.B.31) 

Specifically  what  Is  sought  Is  w(r,<-)  -  w(r)  for  f(r)  of  Interest. 
For  r  large,  from  (II.B.31) 

w(r,c~)  ■  W(r)  •  ’“1/2  vx(x)  .  (II . B. 32) 

NiMerlcal  Integration  by  finite-difference  methods  Is  formidable 
because  the  flow  component  In  the  time-like  direction  u  Is.  In  successively 
thinner  strips  lying  parallel  to  the  boundary,  alternately  In  the  direction 
of  Integration  (stable)  and  opposite  to  It  (unstable).  Though  the  radial 
flow  Is,  on  net.  In  the  direction  of  Integration,  the  Integration  Is 
marginally  stable  and  no  numerical  results  of  real  use  In  the  hurricane 
problem  are  known  to  the  authors. 

At  Carrier's  suggestion,  George  (Carrier,  Hammond,  and  George  1971) 
and  Dergarabedlan  and  Fendell  (1971)  Independently  but  simultaneously 
applied  the  method  of  weighted  residuals  to  the  boundary  value  problem,  and 
found  that,  except  near  the  eyewall  where  the  method  was  Inadequate,  the 
linear  result  given  In  (II.B.32)  sufficed  for  the  nonlinear  problem  as  well. 
Thus,  for  a  form  like  v  -  A(1  -  x/xQ),  for  which  ♦(^.c)  -  0,  w(r,c-~)  - 
-  (A/*0),  a  small  positive  constant  (equivalent  to  about  0.005  mph  down- 
draft  for  physically  Interesting  values). 
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The  reason  weighted-residual  calculations  fall  near  the  axis,  as 
discovered  by  Carrier  (1971a)  by  modified  Oseen  linearization  and  by 
Burggraf,  Stewartson,  and  Belcher  (1971)  by  semi-numerical  analysis  Is 
that  the  structure  of  the  boundary  layer  so  varies  with  radial  distance 
that  adequate  representation  In  terms  of  one  set  of  orthonormal  polynomials 

'  "  ,  Far  fr0m  the  ax1s  °f  sy-etry,  friction  Is  Important  across 

the  entire  layer  of  thickness  0(v/n),/2;  near  the  axis,  friction  Is 

significant  only  Is  a  small  sublayer  near  the  wall  of  thickness  0(r2v/f),/2 
and  the  remainder  of  the  Inflow  layer  of  ON/n)1'2  thickness  Is  Invtscldly  ' 
controlled.  Still,  for  conditions  of  Interest  In  hurricanes,  (II. B  32) 

Is  everywhere  correct  to  within  a  factor  of  two,  and  often  far  better 
Some  results  are  given  In  Figure  (II. B. 5). 

n'B3  and  B  sSppTu  FRICTI0nal  B0UNDAR¥  LA¥ER 

For  the  frictional  boundary  layer  II  and  throughput  supply  I,  one 

takes  the  following  approximations  as  adequate  for  the  quasisteady  mature 
phase: 

1.  The  Prandtl  and  Schmidt  numbers  are  unity; 

2.  the  hydrostatic  approximation  holds; 

3‘  JIle^0U!Idar^  1ayer  approximation  holds  (derivatives  normal 
t1  iob??:dary  excee?  those  tan9ent1  a  1  to  the  bounda?™ 

th«eeS  KetunSr^)!61  *  the 

4.  the  eddy  transfer  Is  adequately  modeled  by  the  laminar 

IlUXh9rwent  :e1ations  for  diffusion  of  massf  momentum 
and  heat  (as  given  by  Flck,  Newton,  and  Fourier,  respectively) 
except  that  the  augmented  kinematic  (eddy)  viscosity  may  vary 
with  radial  position  (but  not  with  ailal  position) ;  a!7  * 

5.  the  mixture  of  dry  air  and  water  vapor  may  be  taken  as  a 

perfect  gas  with  constant  heat  capacity  over  the  range  of 
temperatures  of  Interest  here.  9 

In  view  of  the  limited  understanding  of  quantitative  formulation  of 
cumulus  convection  and  turbulent  transfer,  on  the  cyclone  scale  of  Interest 
here  these  five  approximations  seem  reasonable.  It  can  then  be  shown  that 
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Fig.  1 1 . B . 5a  Nondimensional  results  for  the  frictional  boundary  layer 
obtained  by  the  method  of  weighted  residuals  are  presented  [from  Dergara- 
bedian  and  Fendell  1971jl].  The  divergence  w(x,  c  -*>  °°)  and  the 
volumetric  flux  6  =  - fQ  ( x ,  0  d^  are  presented  for  the  impressed  swirl 

h*  =  1  -  x/x^ ,  x^  =  20,  believed  pertinent  to  a  hurricane  outside  the  eye- 
wall.  Except  near  the  axis  where  nonlinear  inertial  effects  dominate, 
the  linear  Ekman  layer  result,  w(x,  c  -  °°)  -  yx/21/2,  is  an  excellent 
approximation  to  the  numerical  results.  The  volumetric  flux  6 ( x )  is  thus 
linearly  proportional  to  (x^  -  x)  to  good  approximation.  Normalized 
residuals  indicate  large  errors  for  x  <  3,  and  discount  the  premature 
eruption  as  an  artifact  of  the  method,  as  explained  by  Carrier  (l°7la). 

An  improved  solution  for  small  x  indicates  the  adequacy  of  the  linear 
Ekman  solution  to  within  a  factor  of  two.  Since  v  =  (1/75)  mi2/hr  and 
=  (1/16)  rad/hr,  dimensionally  the  results  imply  the  boundary  layer  is 
of  thickness  O(vAa)^2  =  0(1  mi),  the  downflux  into  the  boundary  layer  is 
(2vft)^2  w(x,  ;-+»)=  0(5  x  10  3  mph),  and  the  volumetric  flux  erupting 
up  the  eyewall  is  (2n2  tq2  v/rt)1/2  6(x)  =  0(7  x  103  mi3/hr)  where  v 
characterizes  the  eyewall  angular  momentum  per  unit  mass.  The  implication 
is  that  the  fluid  initially  in  the  boundary  layer  sustains  the  hurricane 
for  about  a  week  and  the  fluid  stored  above  the  boundary  layer  (with 
supplementary  replinishment  of  moisture  from  the  ocean)  can  readily  sustain 
the  hurricane  for  more  than  a  week  more. 
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residual  results  are  erroneous  at  x  <  3,  but  for  larger  x  correctly  indicate 
both  the  reversing  radial-flow  profiles  (with  net  flux  inward)  and  also  the 
swirl  overshooting  its  outer  asymptotes.  The  maximum  magnitude  of  the  radial 
component  is  about  one-third  that  of  the  impressed  swirl  at  a  typical  radial 
position. 


/ 


the  following  equation,  a  generalization  of  ones  given  by  Crocco  In  fluid 

dynamics  and  Shvab  and  Zel'dovlch  In  combustion,  holds  In  the  meteorological 
context  of  Interest  here: 


u 


1H  3H 

3r  w  3z 


a2H 

V  - x  ♦ 


(II .B. 33) 


where  R  denotes  radtatlonal  loss,  taken  by  Carrier,  Hanmond  and  George  to 
be  representable  as 

R  »  -  f(z)  H  , 

with  f(z)  (to  be  discussed  below)  known.  The  definition  of  H  has  been 
given  In  (II. B. 7).  The  boundary-initial  conditions  are  taken  to  be 

z  -  0:  «  ■  Hs;  z  ■  z,:  H  -  H]  ;  (II.B.34*) 

r*  r0!  H  E  Ha<*)  •  (II .B .34b) 


where,  again,  rQ  Is  the  outer  edge  of  the  storm  and  z]  Is  the  top  of  the 
storm.  The  quantities  Ha(z),  H,  -  Ha(z,)  «  Hfl(z  -  0),  and  H  -  H  (z  -  0) 
are  all  given.  Specifically  Hfl(z)  is  the  ambient  distribution  of! 
effectively,  the  equivalent  potential  temperature  times  c  .  Taking 

Hs  =  Ha(z  =  °>  1mP11es  tha*  for  r  >  R  (where  R  Is  the  radius  of  the  eyewall), 
the  ocean  surface  temperature  Is  a  constant  and  the  water  vapor  mass  fraction 
(which  takes  on  Its  saturation  value  at  the  nominally  plane  sea  surface 
z  =  0)  Is  virtually  Independent  of  pressure.  Since  the  thermal  conductivity 
of  water  greatly  exceeds  that  of  air,  uniform  sea-surface  temperature 
appears  to  be  a  good  approximation.  That  H]  =  H  (z  ■  0)  follows  from  the 
definition  of  the  lid  on  the  storm. 

The  terms,  from  left  to  right,  In  (II.B.33)  represent  radial  advectlon- 
axial  convection;  turbulent  diffusion  and  cumulus  convection  (both  para¬ 
meterized  In  v,  which  will  henceforth  be  treated  as  constant,  though  this 
Is  not  necessary);  and  radiation  loss.  Throughout  I,  and  In  II  at  r  =  r 
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since  f(xQ)  ■  0  by  choice  [cf.  (II.B.29)]  and  since  w(3H/3z)  can  be 
neglected  as  explained  below, 


a2H 

f(z)  H  ;  (II.B.35) 

3z‘ 

l.e.,  f(z)  Is  chosen  to  permit  H(z)  ■  H#(z),  given.  Carrier,  Hammond,  and 
George  (1971)  at  this  point  Invoke  Oseen  linearization  arguments.  First, 
they  take  radial  advectlon  as  uniformly  negligible  next  to  axial  convection 
so  the  solution  Is  only  parametrically  dependent  on  r.  This  Implies  slow 
change  of  H  with  r,  and  Is  confirmed  a  posteriori  by  substitution  of  the 
solution  Into  (II.B.33)-(II.B.34).  In  fact.  If  one  takes  H(r,z)  ■  Ha(z), 
which  obeys  the  boundary  and  Initial  conditions,  and  substitutes  this  Into 
the  equation,  one  finds  u(3H/3r)  ■  0;  for  a  typical  value  of  z 

w(aH /az)  , 

— p  -y-  <  0(10  ')  ;  (II  .B. 36) 

v(a  n  /3z  J 


and 


v  — t  -  f(z)  H  (II  .B .37) 

3ZC  # 

by  definition.  The  function  w(r,z)  Is  available  from  results  from  Section 
II. B.l, which  Indicate  that  (w)max  »  w(r,c-*«)  *  W(r),  the  value  at  the  outer- 
edge  of  the  boundary  layer.  In  the  eyewall  w  Is  Increased  by  simple 
continuity  considerations,  to  two  orders  of  magnitude  larger  than  Its 
maximum  value  In  I  or  II  at.  r  >  R.  Hence,  In  the  eyewall,  vertical  con¬ 
vection  dominates, 


w  |y  =  0  ,  (II.B.38) 

and  the  relevant  boundary  data  for  this  hyperbolic  suboperator  Is  that 
given  at  z  =  0  In  (II. B. 34a). 


30 


All  this  discussion,  much  more  carefully  argued  In  Carrier,  Hammond, 

and  George  (1971),  leads  to  two  simply  stated  but  exceedingly  Inportant 
results: 


rQ  >  r  >  R:  H(r,2)  ■  H(rQ,2)  «a(2)  ;  (11.B.39) 

R>r>rfi:  H(r,2)  *  H#(2  ■  0)  ;  (II.B.40) 

where  rQ  Is  the  outer  radial  extent  of  the  eye.  These  last  two  equations 
state  that  for  the  mature  hurricane: 

1.  throughout  the  frictional  boundary  layer  and  the  throughout 
supply,  outside  the  eyewall,  to  within  a  ten  percent  error, 
the  total  stagnation  enthalpy  Is  fixed  at  Its  ambient  ' 
stratification.  Furthermore,  the  ten-percent  correction  Is 
readily  seen  to  be  a  decrease  of  H  with  2  such  that  the 
enthalpy  gradient  at  2  ■  0  Is  Increased  slightly.. 

2.  In  the  eyewall,  the  total  stagnation  enthalpy  Is  constant 
at  Its  sea-level  value,  l.e.,  the  air  Is  rising  on  a  moist 
adlabat. 

Numerical  values  of  Interest  are  the  ambient  net  sea/air  enthalpy 
transfer  and  the  eddy  viscosity  (Carrier,  Hamnond,  and  George  1971): 

-pv(3H/32)|2  ,  0  •  -pv(3Ha/32)|2  .  Q  »  1 .8  X  105  ergs/cm2sec  ; 

(II.B.41) 

v  ■  2.7  X  105  cm2/sec  (II.B.42) 

II. C  THE  RIEHL-MALKUS  MODEL 

The  alternate  theory  of  hurricane  maintenance  Is  that  "...  It  Is 
postulated  that  lowering  of  surface  pressures  In  hurricanes  arises  mainly 
through  an  'extra'  oceanic  heat  source  In  the  storm's  Interior"  (Malkus 
and  Riehl  1960,  p.  12).  To  understand  why  sjch  an  'extra'  oceanic  heat 

source  must  be  postulated  by  Riehl  and  Malkust  one  must  reconstruct  their 
logic. 

First,  the  existence  of  a  frictional  inflow  boundary  layer  is 
acknowledged:  "The  Inflow  Into  a  hurricane  Is  confined  mainly  to  low 
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levels.  Subcloud  air  Is  accelerated  Inward  along  spiral-shaped  trajectorle  ; 
acceleration  results  from  excess  work  done  by  pressure  gradients  over  fric¬ 
tional  retardation"  (Malkus  and  Rlehl  1960,  p.  3).  "...  The  shearing  stress 
vanishes  at  the  top  of  the  Inflow  layer  In  accord  with  the  hypothesis  that 
sv/jz  Is  very  weak  above  the  ground  l^yer"  (Ib1d.,p.  4).  In  a  typical 
hurricane  the  boundary  layer  Is  taken  as  1.1  km  thick,  with  the  eye  boundary 
at  25  km  and  the  extent  of  storm  about  500  km  In  the  Atlantic  and  800-1000  km 
In  the  Pacific.  Malkus  and  Rlehl  choose  to  let  the  relative  vortlclty, 
rather  than  relative  velocity,  vanish  at  the  outer  edge  --  the  result  Is  an 
open  system  with  radial  Inflow  at  the  outer  edge,  as  opposed  to  the  closed 
system  preferred  by  Carrier  and  his  co-workers.  In  a  sample  moderate 
hurricane  Malkus  and  Rlehl  calculate  an  efflux  out  of  the  boundary  layer 
from  the  outer  edge  of  the  eye  out  to  500  km,  an  efflux  that  Increases  with 
decreasing  radius:  "an  average  ascent  rate  of  about  30  cm/sec  or  1  km/hr 
Is  required  at  the  top  of  the  Inflow  layer"  (Ibid.,  p.  7).  This  result 
supposedly  holds  even  though  the  Impressed  swirl  decreases  with  Increasing 
radius;  hence  the  result  contradicts  both  linear  and  nonlinear  theories  for 
the  surface  frictional  layer  under  a  swirling  flow.  There  Is  no  demonstration 
or  Internally  consistent  dynamics  for  the  open  model  of  Malkus  and  Rlehl. 

The  point  of  controversy  Is  that  the  flux  through  the  boundary  layer  does 
not  sink  down  Into  the  boundary  layer  (as  In  the  Carrier  model),  but  rather 
flows  radially  Inward  from  the  outer  edge. 

The  boundary  layer  air,  according  to  Malkus  and  Rlehl,  undergoes 
adiabatic  expansion  as  It  spirals  Inward  toward  the  center,  yet  It  remains 
Isothermal.  This  requires  a  vast,  rather  localized  Increase  In  sea-to-alr 
transfer  between  the  ocean  and  the  contiguous  atmosphere  (In  contrast  to  the 
more  spread  out,  unaugmented  transfer  pictured  In  Carrier's  theory).  That 
the  gradient  normal  to  the  air/sea  Interface  of  temperature  and  of  water 
vapor  mass  fraction  Is  large  enough  to  be  consistent  with  vastly  Increased 
air/sea  transfer  Is  an  article  of  faith: 

In  the  outskirts  of  a  hurricane  the  temperature  of  the  Inflowing 
air  drops  slowly  due  to  adiabatic  expansion  during  (horizontal) 
motion  toward  lower  pressure.  It  Is  one  of  the  remarkable  obser- 

ifnnn°^u  in  hurr*canes  that  thi s  drop  ceases  at  pressures  of  990- 
1000  mb  and  that  thereafter  isothermal  expansion  takes  place. 

Presumably,  the  temperature  difference  between  sea  and  air  attains 
a  value  large  enough  for  the  oceanic  heat  supply  to  take  place  at 
a  sufficient  rate  to  keep  the  temperature  difference  constant. 

(Malkus  and  Riehl  1960,  p.  9). 
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The  actual  transports  [between  sea  and  air],  of  course,  are  very 
W"  ,tf]e  hurricane  compared  to  the  trades.  Sensible  heat  pickup 
is  7Z0  cal/cm  /day,  and  Increase  by  a  factor  of  50  over  the  trades..; 
latent  heat  pickup  Is  2420  cal/cmvday,  higher  by  a  factor  of  12-13 
(Mai kus  and  Rlehl  1960,  p.  12). 

Thus  to  achieve  an  extreme  storm  In  the  framework  of  this  model, 

L turbulent]  transfer  coefficients  enhanced  by  a  factor  of  3-4  aDDear 
to  be  necessary.  (Malkus  and  Rlehl  1960,  p.  16).* 


The  Rlehl-Malkus  theory  that  greatly  augmented  heat  and  mass  transfer 
sustains  the  tropical  cyclone  has.  In  fact,  been  parameterized  Into  all 
existing  computer  simulations  [which  Is  why  all  these  simulations  are  here 
lumped  Into  one  category].  For  example,  It  Is  Interesting  to  note  how 
closely  the  author  of  a  well-known  computer  simulation  (Rosenthal  1971b) 


286-287)'  ar9um*"ts  are  n,ade  at  nore  1e"Sth  elsewhere  by  Rlehl  (1954,  pp. 

Merjy  Published  records,  notably  those  by  Deppemann...  have  proved 
that  the  surface  temperature  outside  the  eye  Is  constant  or  decreases 
very  slightly  toward  the  center.  The  ImplIcatloSs  of  !hls  rmar'Sb'l 

tenrorratur.  ni*.!!11*  noi1ce  “J411  b^ers. .  drew  attention  to  It.  The 

If  2H^h!«“rfaCe  3  r  sp,ral1n9  toward  a  center  should 
decrease  If  adiabatic  expansion  occurred  during  pressure  reduction 

nert  2^  a1  r  enJer1?9  Jhe  circulation  with  the  average  p?o- 
*rKt  eS  ?Lthe.mean  tr°Plcal  atmosphere  should  reach  the  930  mb 

ki  b  iu?1th  3  of  20-5°C  and  specific  humidity  of  17g/ 

kg.  Because  of  condensation,  a  dense  fog  should  prevail  at  the 
ground  Inward  from  the  970  mb  Isobar.  But  this  Is  never  observed 
f°  !°WS  that  -t-he  potential  temperature  of  thp  surf*™  alr  Increases 
^trajectories.  We  also  know  that  the  specific - 1 

hund?edyaJ3  IMo'fee?  th'  Cl°Ud  baSeS  reffla1n  between  a  few 
ul  t^'v^l't^r^yL^Ilf?.  b.°th  1atent  and  ^  Purina. 

1sSQroatlv°Ia1t^Hat^Hdim0lStUre  iIcre;ent  1s  obvious.  The  ocean 
Ilr9rS  fho  and  1arge  amounts  of  water  are  thrown  into  the 

anH  Ik  the.Jorm  °f  spray.  It  Is  hard  to  say  where  the  ocean  ends 
and  where  the  atmosphere  begins!  As  the  air  moves  toward  lower 
pressure  and  begins  to  expand  adlabatlcally,  the  temperature 
ai f ference  between  ocean  and  air  suddenly  Increases.  Since  the 
surface  of  contact  between  air  and  water  Increases  to  many  times 

3rea  °f  Jhe  f*0™’  raP1d  transfer  of  sensible  and 
latent  heat  from  ocean  to  air  Is  made  possible.  In  the  outskirts 

say  beyond  the  990  mb  Isobar,  the  tumoll  is  less and  the ^  process 
of  heat  transfer  Is  not  operative.  (Riehl  1954,  pp.  286-287). 
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reflects  the  Riehl-Malkus  theory  In  his  most  recent  publication  and  how  he 
notes  similar  logic  in  the  work  of  another  computer  modeler  of  tropical 
cyclones  (Ooyama  1969): 


{ 


Air-sea  exchanges  of  sensible  and  latent  heat  have  long  been 
considered  important  Ingredients  In  the  development  and  main¬ 
tenance  of  tropical  storms.  Palmen  (1948)  showed,  on  a 
climatological  basis,  that  tropical  storms  form  primarily  over 
warm  ocean  waters  (Teea  >  26°C).  Malkus  and  Rlehl  (I960)  showed 
that  the  deep  central  pressures  associated  with  hurricanes  could 
not  be  explained  hydrostatlcal ly  unless  the  equivalent  potential 
temperature,  e-,  in  the  boundary  layer  was  10°  to  15°K  greater 
than  that  of  the  mean  tropical  atmosphere.  Byers  (1944)  pointed 
out  that  the  observed  near-isothermal  conditions  for  inward 


spiraling  air  in  the  hurricane  boundary  layer  required  a  source 
of  sensible  heat  to  compensate  for  the  cooling  due  to  adiabatic 
expansion. . . 


Ooyama  (1969)  found  drastic  reductions  In  the  strength  of  his 
model  storm  when  the  air-sea  exchanges  of  sensible  and  latent 
heat  were  suppressed.  He  pointed  out  that  at  sufficiently  large 
radii,  the  boundary  layer  is  divergent  (the  so-called  Ekman 
layer  "sucking")...  This  subsidence  tends  to  decrease  the  boundary 
layer  e„  since  ae  /3z  <  0  In  the  lower  troposphere.  Ooyama 
argued  that  unless  the  energy  supply  from  the  ocean  can  again 
raise  the  ee  of  the  boundary  layer  air  to  sufficiently  large 
values  before  the  Inflowing  air  reaches  the  Inner  region,  and 
the  convective  activity  will  diminish  In  those  regions  and, 
hence,  the  storm  will  begin  to  weaken. 

Ooyama 's  line  of  reasoning  can  be  extended  to  show  that  evaporation 
is  far  more  important  than  sensible  heat  flux.  The  air  sucked  into 
the  boundary  layer  has  a  higher  potential  temperature  than  the 
original  boundary  layer  air.  The  subsiding  air  has  a  smaller  ee 
only  because  it  is  relatively  dry.  (Rosenthal  1971b,  p.  772). 


An  extensive  presentation  has  been  given  to  the  Riehl-Malkus  theory, 
and  an  extensive  criticism  will  now  follow,  because  the  theory  has  such 
wide  acceptance.  For  example,  a  recent  National  Science  Foundation  document, 
in  reviewing  the  work  of  Ooyama  at  NYU  —  work  emphasizing  the  reliance  on 
high  sea  temperature  and  attendant  large  enthalpy  transfer  from  sea  to 
air  --,  states:  "Besides  the  fact  that  the  model  has  succeeded  in  simu¬ 
lating  many  important  aspects  of  the  hurricane,  it  has  also  demonstrated 
the  importance  of  oceanic  latent-heat  supply  to  both  the  development  and 
maintenance  of  the  tropical  storm"  (NSF  1969,  p.  96).  Battan  (1969,  p.  117) 
cites  Ooyama 's  results  with  no  reservations:  "Ooyama1 s  research  shows  that 
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storm  dynamics  depends  mostly  on  the  temperature  of  the  ocean  water  under 
the  hurricane."  In  fact,  the  very  title  of  a  recent  Russian  article 
citing  Riehl  and  Malkus  papers  suggests  widespread  acceptance:  "The 
Power  of  a  Tropical  Cyclone  as  a  Function  of  the  Underlying  Sea  Surface 
Temperature"  (Shuleykin  1971).  The  fact  is  that  ocean  temperatures  are 
difficult  to  determine  from  currently  available  records,  and  that  hurri¬ 
canes  cause  upwelling  of  lower,  colder  water  to  confuse  matters  further 
(Perlroth  1967).  Gentry  (1969,  p.  406)  presents  data  relating  "...  the 
maximum  intensity  of  several  tropical  cyclones  to  the  temperatures  of  the 
sea  beneath  them  and  shows  that  both  severe  and  weak  tropical  cyclones 
occur  when  ocean  temperatures  are  relatively  high.  This  suggests  that 
variations  in  parameters  other  than  the  transfer  of  heat  from  the  ocean  to 
the  atmosphere  also  influence  the  storm's  intensity...".  Actually  Gentry's 
data  [Figure  (II.B.6)]  indicate  that  many  intense  tropical  cyclones  He 
over  relatively  cold  ocean  water  {<  28.0°C).  Brand  (1971)  cites  a  supposed 
correlation  of  central  pressure  deficit  with  sea-surface  temperature 
throughout  the  lifespan  of  Hurricane  Esther  of  9-26  September  1961.  But 
during  one  twelve-hour  period  of  constant  sea  temperature,  the  central 
pressure  rose  15  mb;  and  for  three  days  while  the  sea  temperature  hovered 
about  84°F,  the  central  pressure  nonmonotonically  rose  from  930  to  955  mb; 
and  when  the  sea  temperature  was  at  86°F,  at  various  times  the  central 
pressure  was  as  low  as  927  mb  and  as  high  as  953  mb.  Perlroth  (1967)  shows 
that  for  Hurricane  Ginny  of  1963,  for  five  days  while  the  sea  temperature 
hovered  near  80°F,  the  central  pressure  nonmonotonically  fell  from  995  to 
970  mb  [Figure  II.B.6].  While  some  correlation  of  central  pressure  deficit 
and  sea  surface  temperature  definitely  exists,  tropical  cyclone  intensity 
cannot  be  completely  correlated  In  so  simple  a  fashion.  Thus,  because  of 
the  fact  that  It  Is  a  starting  point  for  almost  all  subseouent  theoretical 
research  on  tropical  cyclones,  there  seems  reason  to  scrutinize  the  Riehl - 
Malkus  theory  carefully. 

There  would  seem  to  be  several  errors  In  the  Malkus-Riehl  reasoning. 
First,  even  for  the  linearized  Ekman  layer  at  the  outer  edge  of  the  tropical 
cyclone.^1 t  Is  Inappropriate  to  use  the  adiabatic-expansion  relation 
P  %  TY  Y  *  which  1mPlies  f°r  steady  motion  that  the  entropy  is  constant 


35 


Fig.  II. B. 6  Above,  the  minimum  sea-level  pressure  measured  for  a 
tropical  cyclone  is  plotted  against  the  local  sea  temperature  [from  Gentry 
(1969a,  p.  406)].  Below,  the  central  sea-level  pressure  is  plotted  as  a  function 
of  the  local  sea  surface  temperature  for  hurricane  Ginn y,  21-28  October  1963, 
which  passed  over  the  Gulf  Stream  on  24  October  [from  Perlroth  (1967,  p.  266)]. 
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R<eh1  ^  P’  2861  “"*«*  ■»«**  that  neglect 
L  1  “  6"C  dr°P  (wUh  condensation)  end  e  7.3°C  drop  (dry) 

fro™  the  mean  autumnal  see-level  ambient  reported  by  Jordan  (1957)  for  the 

West  Indies  area.  In  dropping  from  1014  mb  to  930  mb.  The  Isentroplc 
assumption  -.s  unjustified  and  the  temperature  decrease  thereby  anticipated 


What  temperature  decrease  Is  to  be  expected?  It  has  been  discussed 
earlier  In  Section  II. 8. 3  that  an  approximate  Integral  for  the  surface 
r  ct  onal  layer,  more  than  adequate  for  most  purposes,  Is  given  by 


H(r.t)  =CpT  +  ga  +  LY  +  (q2/2)  J  H(r0.z);  (II.C.l) 

l.e..  the  total  stagnation  enthalpy  H  throughout  the  boundary  layer  between 
the  eyewall  and  the  outer  edge.  Is  approximately  fixed  at  the  stratification 
of  the  ambient  environment  In  which  the  tropical  cyclone  was  generated. 
Again,  ro  denotes  the  outer  edge  of  the  storm.  At  z  =  0.  Y  Is  close  to  Its 
saturated  value  and  this  Is  primarily  a  function  of  the  sea  surface 
temperature,  which  Is  adequately  modeled  as  a  constant.  Hence, 


"o  <r.°+>  ?  cp 


T  + 


q2/2 


(ro.0+) 


const. 


If  q  a  o,  T  -  299.4  at  r 
Is  a  decrement  of  4.0°C. 


V  z  a  0  »  then  at  q  =  200  mph,  T  =  295.4°K;  this 


It  will  now  be  suggested  that  this  Is  the  decrement  Rtehl  should  have 
anticipated.  The  cydostrophlc  balance  [Fletcher  (1955)  formula]  gives 


(II. C. 2) 


f  the  radial  profile  of  the  swirl  above  the  boundary  layer  is  given  by 

(II.B.ll).  Here  n  =  0.6  and  cs  .  1.2  X  10-3  g/cm3  are  adopted.  The 

ambient  sea-level  pressure  p$  Is  taken  as  1014  mb  according  to  Jordan, 
and  Rlehl  for  his  calculations  adopted  a  case  in  which  the  central  sea- 

level  pressure  fell  to  930  mb.  Substitution  In  (II. B.12)  gives  (v)  1 

209  mph,  about  as  large  as  any  value  reliably  reported  for  a  troplcarcydone. 
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Thus,  for  the  pressure  difference  adopted  by  Riehl,  equivalent  to 
(v)max  a  209  mph,  the  appropriate  generalized  Crocco-Shvab-Zel 'dovich-li ke 
integral  gives  a  temperature  drop  of  about  4.3°C;  Riehl's  use  of  the 
adiabatic  relation  led  him  to  expect  a  temperature  drop  at  least  three-halves 
as  large*. 

Riehl  (1954),  Malkus  and  Riehl  (1960),  and  Rosenthal  (1971b)  all  cite 
measurements  that  purportedly  show  that  the  sea-level  atmospheric  tempera¬ 
ture  in  a  mature  tropical  cyclone  is  constant.  "This  remarkable  fact..." 
(Riehl  1954,  p.  286)  Is  based  on  measurements  cited  by  Riehl,  measurements 
made  over  twenty-five  years  ago  when  hurricane  speeds  were  not  commonly 
believed  to  attain  200  mph  levels,  and  when  the  distinction  between  static 
enthalpy  Cp  T  and  stagnation  enthalpy  (Cp  T  +  q^/2)  [normally  unimportant 
in  meteorology  but  often  significant  in  modern  high-speed  aerodynamics] 
was  not  always  carefully  observed.  The  fact  is  that  accurate  sea-level 
measurements  In  Intense  hurricanes  are  not  a  simple  matter  to  this  day 
because  of  reliance  on  make-shift  combinations  of  obsolete  military  air¬ 
craft  and  radar  (Meyer  1971).  The  belief  here  is  that  the  static  tempera¬ 
ture  actually  decreases  as  predicted  by  the  Carrier  theory  when  properly 
measured. + 


*It  is  interesting  to  note  that  according  to  Riehl  (1954,  p.  286)  a  dry- 
adiabatic  expansion  from  1014  to  970  mb  would  cause  condensation.  Equi¬ 
valently,  he  is  stating  that  about  a  3.7°C  drop  would  cause  condensation 
in  sea-level  tropical  air  in  the  hurricane  season.  Since  even  an  intense 
(200  mph)  hurricane  really  gives  a  total  drop  of  only  4.0°C,  one  sees  why 
the  boundary  layer  is  quite  cloud-free  in  as  far  as  the  eyewall. 

+Palmen  and  Newton  (1969,  p.  478)  still  cite,  second-hand,  measurements  over 
thirty- five  years  old  asserting  the  constancy  of  the  temperature  in  the 
surface  frictional  layer.  One  must  recall  that  aircraft  penetration  of 
hurricanes  is  less  than  twenty  years  old.  They  also  cite  a  1954  work  by 
Arakawa  in  which  the  wet-bulb  temperature  held  nearly  constant  during  the 
passage  of  an  898-mb  typhoon  over  a  Japanese  naval  fleet;  this  is  inter¬ 
preted  to  prove  that  the  equivalent  potential  temperature  of  the  surface 
air  rose  25°K  in  spiraling  in  from  the  outskirts  to  the  center,  and  hence 
that  an  internal  heat  enthalpy  source  in  the  form  of  augmented  sea-air 
latent  and  sensible  heat  transfer  is  operative.  However,  reference  to 
even  a  recent  book  on  meteorology  [e.g.,  Hess  (1959,  p.  61)]  will  show  that 
the  dynamic  contribution  q2/2  to  the  total  enthalpy  of  the  gas  is  neglected 
in  computing  the  temperature  from  the  measured  wet-bulb  temperature;  this 
normally  insignificant  correction  makes  only  a  little  over  a  one-percent 
contribution  to  the  total  enthalpy  even  in  a  hurricane.  But  confusion  about 
this  small  effect  has  had  a  profound  influence  on  the  evolution  of  tropical 
cyclone  theory.  Actually,  radially  constant  wet-bulb  temperature  confirms  (II.C.l). 
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Thus  one  is  led  to  the  conclusion  that  Riehl  believed  (1)  the 
hurricane  draws  upon  local  ambient  air  (not  warm  moist  air  sinking  into 
the  boundary  layer);  (2)  a  large  thermal  fall  occurs  while  spiraling 
through  the  boundary  layer  (at  least  50%  above  what  seems  correct);  and  (3) 
the  isothermal  flow,  purportedly  measured  in  the  boundary  layer,  occurs 
(whereas  it  almost  surely  does  not).  The  result  is  that  Riehl  decided 
that  an  internal  oceanic  heating  source  was  necessary  and  postulated  greatly 
augmented  heat  and  mass  transfer  from  ocean  to  atmosphere  in  a  hurricane, 
especially  as  the  center  is  approached  (whereas  the  result  compatible  with 
the  conservation  laws  is  that  the  increase  in  latent  and  sensible  heat 
transfer  where  a  hurricane  lies  over  what  would  be  transferred  in  its 
absence  is  at  most  five  or  ten  percent,  and  this  slight  augmentation  occurs 
at  the  outer  edge  and  the  slight  augmentation  falls  off  as  the  center  is 
approached).  This  basic  slip  is  incorporated  into  every  numerical  model 
known  to  the  authors  and  raises  questions  about  many  conclusions  drawn  from 
the  extensive  programming,  computer  solution,  and  interpretation  that  has 
gone  on  for  several  years  now. 

Rosenthal  (1971b),  in  response  to  publications  of  the  authors,  now 
concurs  that  greatly  augmented  sensible  heat  transfer  is  dispensible,  but 
remains  resolute  that  greatly  augmented  latent  heat  transfer  from  the  ocean 
to  the  atmosphere  -  critical  to  his  theory  and  to  Ooyama's  --  is  necessary 
for  hurricane  generation  and  maintenance.  His  statements  that  higher-level 


The  assertion  by  Riehl  (1954,  p.  287)  that  much  spray  is  tossed  into  the 
air  within  hurricanes  does  not  in  itself  assure  augmented  net  enthalpy 
transfer  from  sea  to  air,  since  heat  must  be  drawn  from  the  air  to  evaporate 
the  drops  for  later  condensation  of  the  water  vapor  in  the  eyewall  and  inner 
rainbands .  Incidentally,  one  must  also  be  careful  about  extrapolating 
empirical  laws  relating  ocean  evaporation  to  wind  speed  to  near-saturation 
high-speed  conditions  for  which  they  were  not  devised.  Finally,  the  greatly 
aUjniMn^f^  anc*  sei^s i bl e  heat  sea-to-air  transfer  postulated  by  Riehl 

and  Malkus  (I960,  p.  17)  amounts  to  but  about  three  percent  of  the  total 
enthalpy  of  the  inflowing  tropical  air.  Yet,  even  acknowledging  this,  they 
insisted  that. this  marginal  amount  is  critical  to  both  genesis  and  mainte- 
nance  of  hurricanes.  The  Carrier  model  states  that  the  increase  postulated 
by  Riehl  and  Malkus  is  an  order  of  magnitude  too  large  and  is  probably  not 
the  important  critical  factor.  J 
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air  sinking  into  the  boundary  layer  is  relatively  dry  and  of  lower 
equivalent  potential  temperature  may  arise  from  misunderstanding  about  the 
quasisteady  nature  of  the  Carrier  model  and  supporting  analysis  for  the 
mature  hurricane.  The  air  slowly  sinking  into  the  boundary  layer  at  about 
0.005  mph  is  enriched  by  plume  and  radiative  transfer  of  water  vapor  and 
heat  from  the  ocean  to  atmosphere.  These  ambient  mechanisms,  as  has  been 
repeatedly  emphasized,  persist,  neither  augmented  nor  decremented  to  any 
important  degree,  within  the  hurricane;  they  help  to  compensate  for  rain- 
out  in  the  outer  spiral  bands.  If  Indeed  the  ambient  tropical  profile  for 
equivalent  potential  temperature  does  persist  with  little  change  in  the 
hurricane,  much  of  the  sea/air  transfer  continues  to  pass  across  the 
boundary  layer  with  little  dimunition,  just  as  in  the  ambient;  the  enthalpy 
enriches  the  air  in  the  700-900  mb  strata  such  that  by  the  time  this  air 
sinks  into  the  boundary  layer,  it  is  much  like  the  air  originally  in  the 
900-1000  mb  strata.  Only  as  the  hurricane  leaves  the  tropical  oceans  is 
this  normal  transfer  reduced;  the  air  entering  the  boundary  layer  eventually 
is  of  lower  equivalent  potential  temperature  since  it  comes  from  air  which 
is  originally  higher  in  the  tropical  ambient  (hence  colder  and  drier)  and 
which  is  not  appreciably  enriched  as  it  descends.  In  this  way  traverse 
over  ocean  patches  of  varying  temperature  can  cause  the  well-known  non¬ 
monotonic  perturbations  in  hurricane  intensity  within  the  general  level  of 
strength  computed  from  the  spawning  ambient  as  discussed  earlier. 

If  the  basic  overall  physics  of  all  computer  models  is  the  same,  why  is 
there  a  proliferation  of  programs  (Ooyama  1969;  Rosenthal  1970;  Sundqvist 
1970;  Yamasaki  1968;  Kurihara  1971)?  The  reason  is  that  while  the  programs 
all  agree  with  Riehl-Malkus  concepts  on  the  cyclone  scale,  there  are 
physical  processes  that  are  difficult  to  parameterize  (turbulent  diffusion, 
radiational  transfer);  also,  there  is  the  question  of  how  to  model  the 
cumulus-convection  scale  within  the  cyclone-scale  program. 

First,  it  is  questionable  whether  it  is  feasible  to  seek  a  uniformly 
valid  solution  to  the  entire  cyclone  when  large  gradients  occur  over 
relatively  small  scales  in  important  subregions  of  the  storm  (e.g.,  the 
frictional  boundary  layer  and  the  eyewall),  while  small  gradients  occur 
over  relatively  large  scales  in  the  bulk  of  the  storm  (the  rapidly  swirling 
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regions  and  the  outflow  layer  aloft).  On  current  computers  usually  a  more- 
or-less  fixed  grid  of  10  km  or  20  km  radial  resolution  and  at  most  thirteen 
layers  of  vertical  resolution  is  adopted  by  practical  consl derations;  the 
overwhelming  bulk  of  the  grid  points  then  lie  outside  the  eye,  the  eyewall, 
and  the  frictional  boundary  layer  —  where  Important  processes  are  occurring.* 

The  domain  sizes  are  quite  limited  {typically  440  km)  so  relatively  weak 
boundary  conditions  (requiring  only  that  purely  advective  influx  occur  at 
the  side  boundaries)  are  employed  arid  much  of  the  storm  lies  outside  the 
domain  of  computation.  When  a  closed  system  Is  studied  (no  radial  inflow 
over  side  boundaries),  the  peak  Intensity  increases  monotonically  with 
domain  size,  even  to  1200  km  (Rosenthal  1971b).  Oust  what  can  be  discerned 
from  such  results  Is  mute.  It  seems  more  useful  to  subdivide  the  storm 
Into  natural  portions  where  different  gradients  and  phenomena  are  operative, 
as  In  Carrier's  approach. 

Next,  what  may  be  proven  with  current  numerical  models  deserves 
consideration.  The  air/sea  interaction  postulated  by  Riehl  and  Malkus  has 
been  parameterized  Into  Rosenthal's  model  and  into  Ooyama's,  in  lieu  of  a 
solution  of  the  boundary  layer.  Naturally  results  from  both  models  reflect 
this  formulation  and  do  not  corroborate  its  physical  validity,  although  Rosenthal 
would  disagree  (Rosenthal  1971b,  p.  767,  771-772).  In  fact,  the 
computation  Is  curve-fitting  In  the  sense  that  numerical  techniques  are 
rated  "...  on  intuitive  meteorological  inspection  of  test  results"  (Anthes, 
Rosenthal,  and  Trout  1971,  p.  747).  Errors  introduced  by  finite-differencing 
are  used  to  simulate  physical  phenomena  (as  lateral  mixing);  the  more 
accurate  the  differencing,  the  worse  the  results.  Deterioration  of  results 

♦Even  devoting  four  points  to  defining  a  frictional  boundary  layer  profile 
Is  quite  marginal . 

+"The  deterioration  of  the  solutions  with  the  introduction  of  the  centered 
difference  scheme  was  not  anticipated,  and,  indeed,  was  quite  disappointing. 

Not  only  does  the  less  accurate  upstream  method  provide  model  storms  with 
more  acceptable  structure  and  better  consistency  between  wind  and  pressure 
but  also  (sic)  consistency  provides  an  Internal  dissipation  of  kinetic 
energy  of  the  same  order  of  magnitude  as  the  suiiace  dissipation...  It 
appears...  that  this  is  the  correct  proportionality  between  internal  dissi¬ 
pation  and  surface  dissipation...  These  'beneficial'  aspects  of  upstream 
differencing  are  clearly  fortuitous...  They  seem  to  lead  to  the  conclusion 
that,  with  our  present  lack  of  knowledge  concerning  the  interactions  (continued) 
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with  increased  precision  is  common  In  curve-fitting.  It  also  appears 
that  numerical  models  are  so  sensitive  to  small  changes  in  Initial  data 
and  in  parameterizatlons  of  frictional  effects  (necessitated  by  the  absence 
of  a  solution  for  the  boundary  layer)  that  little  can  be  learned  about 
intensification  and,  at  most,  only  the  mature-stage  structure  is  meaningful.* 

II. D  IMPLICATIONS  OF  THE  MODELS  ON  SEEDING 

The  current  mode  of  seeding  is  to  introduce  silver  iodide  crystals  In 
supercooled  water  believed  to  exist  high  In  the  eyewall,  and  sometimes  also 
in  rainbands  both  close  to  and  far  from  the  eyewall.  The  silver  iodide 
will  hopefully  cause  the  water  droplets  to  freeze  and  release  the  heat  of 
fusion  (Battan  1969).  The  proposed  mechanism  by  which  this  heat  release 
{and  attendant  temperature  rise  and  density  decrease)  causes  an  amelioration 
of  tropical  cyclone  intensity  has  been  altered  several  times  (Rosenthal 
1971a;  Gentry  1971b)  and  apparently  is  now  uncertain  [Gentry  (1971a)  proposes 
several  mechanisms];  even  in  its  original  concept  (Simpson  and  Malkus  1964), 
the  mechanism  by  which  seeding  was  to  be  efficacious  seemed  nebulous 
to  the  current  authors.  The  one  accepted  point  is  that  seeding  in  the 
nascent  eye  of  a  developing  tropical  storm  should  be  avoided  since  this 
procedure  would  probably  abet  intensification  (Rosenthal  1971a). 


h 

(continued)  between  the  cumulus  scale  and  the  macroscale,  the  diffusive 
effects  provided  by  upstream  differencing  are  probably  as  good  a  representa¬ 
tion  of  the  statistical  effect  of  the  cumulus  motions  on  the  macroscale 
velocity  fields  as  anything  currently  available.  Such  a  conclusion,  of 
course,  only  points  to  a  high  degree  of  ignorance  with  regard  to  an 
extremely  important  meteorological  problem.  It  is  by  no  means  a  solution." 
(Rosenthal  1970,  p.  657-658). 

*"7he  'organizational'  period  is  about  twice  as  long  as  that  found  in  our 
previously  published  results...  This  is  primarily  a  result  of  replacing 
the  constant  drag  coefficient  (3  X  10~3)  with  the  variable  Cp... 

The  time  needed  for  the  model  cyclone  to  become  organized  is  also  highly 
sensitive  to  the  arbitrary  initial  conditions... 

The  material  presented  in  the  last  few  paragraphs  indicates  that  the  length 
of  the  organizational  period,  as  given  by  model  calculations,  is  only  of 

significance  when  experiments  are  compared  against  each  other."  (Rosenthal 
1971 ,  p.  769. ) 
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S-'nce  the  storm  is  naturally  oscillating  in  intensity  and  the  threat 
of  litigation  has  constrained  the  number  of  seeding  experiments,  there  is 
little  to  permit  discrimination  of  natural  and  artificially  induced  changes 
in  intensity.  This  is  aggravated  by  the  fact  that  even  proponents  of 
seeding  anticipate  only  a  ten- to-fif teen  percent  decrease  in  maximum  winds. 
In  the  one  case  in  which  larger  decreases  were  noted,  the  anomaly  is  now 
attributed  to  synoptic  peculiarities  relating  to  the  upper-level  outflow 
(Hawkins  1971). 

The  National  Hurricane  Research  Laboratory  (Gentry  1969a;  Gentry  1969b; 
Gentry  1970)  has  discussed  a  six-to-twel ve  hour  cycle  of  amelioration  after 
seeding;  physical  basis  for  this  time  scale  has  yet  t j  come  forth.  Further, 
if  the  central  pressure  deficit  is  reduced  as  reported,  eyewall  seeding 
must  alter  the  eye  in  an  as-yet  unidentified  manner.  More  complete  post- 
seeding  probing  of  the  tropical  cyclone  would  be  helpful  in  evaluating 
these  claims. 

Rather  similar  computer  models  have  produced  different  guidance  with 
regard  to  current  seeding  practice.  Sundqvist  (1971)  states  flatly  that 
it  will  intensify  the  hurricane.  Rosenthal  (1971a)  suggests  that  it  will 
displace  the  maximum  winds  to  greater  radii,  reduce  the  maximum  winds  by 
ten- to-fif teen  percent*,  and  increase  outer  winds  by  ten  to  fifteen  percent. 
The  reasons  for  the  discrepancy  are  not  available  because  some  details 
remain  unpublished. 

The  present  authors  have  been,  and  still  are,  skeptical  about  the 
effectiveness  of  current  seeding  practices.  If  the  Carrier  model  is  valid, 
silver  iodide  seeding  can  only  transiently  upset  the  stable  hurricane 
configuration.  Under  this  model,  warm-fog  dispersal  methods  would  have  to 
be  applied  to  the  entire  "throughput  supply"  layer  of  warm  moist  air  to 
achieve  the  significant  goal  of  premature  rainout  in  outer  spiral  bands. 

The  layer  of  warm  moist  air  is  so  spatially  extensive  that  such  attempts 
seem  somewhat  impractical.  Seeding  to  divert  the  path  slightly  holds 


♦Statistical  treatment  of  a  model  suggests  multimillion  dollar  annual  savings 
in  damage  from  seedings  which  would  reduce  peak  winds  by  fifteen  percent 
(Boyd,  Howard,  Matheson,  and  North  1971). 
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little  better  promise  than  seeding  to  alleviate  intensity,  since  there 
seems  no  way  to  discern  what  path  alterations  were  due  to  human  interven¬ 
tion  under  current  understanding. 
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III.  THEORY  OF  TROPICAL  CYCLONE  INTENSIFICATION 


A  more  difficult  problem  than  understanding  the  quasisteady  structure 
of  a  mature  tropical  cyclone  is  developing  a  transient  analysis  explaining 
intensification  from  a  weaker  disturbance.  The  subtleties  of  the  problem 
are  suggested  by  the  fact  that  of  several  hundreds  of  disturbances  over 
the  tropical  oceans  in  autumn,  annually  but  about  fifty  cyclones  form  and 
there  is  currently  no  way  to  predict  which  disturbances  will  intensify. 

No  complete  picture  of  intensification  will  be  given  here.  But  for 
the  Carrier  model,  the  only  one  deemed  worthy  of  further  study,  the  out¬ 
line  of  an  intensification  theory  will  be  presented,  together  with  several 
successively  more  sophisticated  models  aimed  at  quantitative  exploration 
of  details. 

III. A  CARRIER'S  OUTLINE  OF  INTENSIFICATION 

Carrier  (1571b)  has  already  sketched  an  intensification  process  by 
which  his  quasisteady  mature  model  evolves  in  time  from  a  tropical 
depression.  Tracing  back  to  even  earlier  evolution  seems  premature  at 
the  current  state  of  understanding.  It  can  hardly  be  overemphasized  that 
everything  which  follows  is  either  a  reproduction  or  refinement  of  Carrier's 
earlier  work. 

In  a  tropical  depression  there  is.  radial  inflow  in  region  I  and  there 
is  weak  Rankine-vortex-like  swirling,  tne  maximum  azimuthal  speed  lying 
relatively  far  from  the  axis  (axisymmetric  model).  The  swirling  would 
quickly  establish  a  shear  layer  beneath  it;  there  is  weak  upflux  out  of 
(downflux  into)  the  surface  layer  where  the  swirling  speed  increases 
(decreases)  with  radial  distance.  Since  linear  theory  correctly  predicts 
the  downflux  in  the  mature  stage  with  rapid  swirl,  linear  theory  certainly 
suffices  during  intensification.  However,  since  equilibration  of  the 
boundary  layer  requires  times  of  0(n-1),  or  roughly  sixteen  hours  in  the 
tropics,  a  transient  linear  theory  will  be  required.  It  must  yield  a  radial 
influx  through  the  boundary  layer  II  in  excess  of  the  radial  inflow  speed 
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EQUIVALENT  POTENTIAL  TEMPERATURE,  8g 


THF  DISTRIBUTION  OF  0  AS  A  FUNCTION  OF  PRESSURE  (p)  JN  MIUI- 
I”L  nJ  HFirHT 7,1 IN  *  KILOMETERS  FOR  DAYS  WITH  FOUR  CHARAC¬ 
TER  is  nr  WFATHER  TYPES-  AVERAGE  (SOLID);  SUPPRESSED  CONVECTION 
DOT  ID) MoS  RA ELY  ENHANCED  CONVECTION  (LONG  DASHES);  AND 

$1  WDNGLY^^NH^WC  ED  CONVECTION  |SHO*rD*SHeSl.JH£OIMMtNC£S 

IN  THE  DISTRIBUTIONS  RESULT  PRIMARILY  FROM  THE  VERTICAL  TRANS 
PORT  OF  WATER  VAPOR  BY  THE  CONVECTIVE  CLOUDS. 


TOTAL  ENTHALPY  H  =  cpT  +  q  ♦  LY  +  g* 

VERSUS  ALTITUDE  i,  WHERE  q  IS  WIND 
SPEED,  L  IS  LATENT  HEAT  OF  VAPORIZATION 
PER  UNIT  MASS  AND  Y  IS  THE  MASS  FRACTION 
OF  WATER  VAPOR.  THE  CURVE  LABELLED 
a  IS  TYPICAL  OF  THE  AMBIENT  ATMOSPHERE 
AND  m  IS  THE  SEQUENCE  OF  STATES  ACHIEVED 
BY  STARTING  AT  SEA  LEVEL  AND  FOLLOWING 
&  MOIST  ADIABATIC  PROCESS  (AS  IN  THE  EYEWALL) 


Fig.  III.A.1  On  the  left  is  a  plot  of  the  equivalent  potential 
temperature  measured  near  Barbados  (in  the  Lesser  Antilles)  in  July-August 
1968,  showing  a  minimum  at  about  700  mb,  which  becomes  less  clearly  defined 
as  convection  increases  [taken  with  caption  from  Garstang,  La  Seur,  Marsh, 
Hadlock,  and  Petersen  (1970,  p.  494)].  On  the  right  are  the  approximations 
used  here;  the  total  enthalpy  H,  since  the  dynamic  contribution  is 
nealiaible  in  the  ambient  and  never  a  dominant  contribution,  is  effectively 
c  o  .  Clearly  the  role  of  enhanced  convection  is  to  alter 
ambient  profile  for  H  in  the  direction  of  a  vertical  line. 


the  same  as 
the  typical 
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Fig.  I II. A. 2  In  this  schematic  view  of  the  flow  configuration  and 
circumferential  velocity  distribution  in  an  intensifying  tropical  depression 
at  some  early  time  t  =  t-j  (say),  the  interface  C-C  between  the  new  and 
initial  air  in  the  core  is  idealized  as  horizontal  for  convenience. 
Meticulous  details  such  as  this  are  of  no  current  concern  since  many  points 
have  yet  to  be  quantitatively  resolved. 


z 


Fig.  III. A. 3  Intensification  from  tropical  depression  to  hurricane 
has  proceeded  to  a  more  advanced  stage  in  this  schematic  diagram,  holding 
at  t  =  t^  >  t-j .  The  magnitude  of  the  maximum  swirl  is  increased  and  its 
position  lies  closer  to  the  axis  of  symmetry. 
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in  I.  Incidentally,  the  times  characterizing  significant  change  in  the 
thermodynamic  state  of  the  lower  tropical  ambient  are  probably  on  the  order 
of  100  hours,  a  time  span  partly  related  to  the  eddy  viscosity  appropriate 
for  the  tropical  atmosphere.  Hence,  unless  the  intensification  time  from 
depression  to  cyclone  greatly  exceeds  100  hours  --  and  this  seems  dubious  -- 
the  thermodynamic  state  of  the  ambient  may  be  considered  fixed  throughout 
intensi fication. 

The  air  erupting  from  the  boundary  layer  begins  to  displace  the  air 
initially  in  the  central  core  of  the  developing  storm.  The  air  initially 
present  in  the  core  is  of  slightly  lower  pressure  than  the  ambient  air  at 
the  edge  of  the  storm,  but  not  vastly  different  in  vertical  stratification. 
The  air  erupting  from  the  boundary  layer  under  the  Rankine- vortex-1  ike 
swirl  displaces  the  air  initially  present  in  the  core  vertically  upward; 
since  there  is  a  "lid"  on  top  of  the  storm,  the  vertically  displaced 
initial  air  is,  near  the  top  of  the  core,  squeezed  radially  outward. 

The  following  competition  develops.  The  new  air  rising  out  of  the 
boundary  layer  is  drawn  entirely  from  relatively  warm  moist  air  near  the 
bottom  of  the  atmosphere.  Thus,  displacing  the  air  initially  present  in 
the  core  is  relatively  light  air.  On  the  other  hand,  the  convective  motion 
of  new  air  is  small,  especially  at  early  times,  and  the  ambient  processes 
(turbulent  diffusion,  radiational  cooling,  cumulus  convection)  try  to 
maintain  the  original,  near-ambient  stratification  in  the  core.  If  the 
convective  displacement  wins  out,  then  the  core  becomes  lighter  and  lighter, 
relative  to  a  column  of  air  at  the  outer  edge  of  the  storm  [see  Figure 
(III.A.l)]. 

Thus,  the  swirling  in  I  has  led  to  a  downflux  in  II,  a  spiraling 
inward  in  the  boundary  layer  and  an  upflux  into  the  core,  and  a  lightening 
of  the  core  by  hydrostatic  considerations.  For  dynamic  consistency,  the 
centrifugal  force  (anticipated  to  be  the  dominant  inertial  effect)  must 
increase  to  balance  the  augmented  radial  pressure  gradient.  Since  angular 
momentum  is  conserved  in  I,  where  friction  is  negligible,  the  fluid  particles 
must  necessarily  move  in  closer  to  the  axis  of  symmetry  (axis  of  rotation). 
The  result  is  that  in  time,  in  the  Rankine-vortex-like  swirl  distribution, 
the  maximum  azimuthal  speed  increases  in  magnitude  and  the  position  of  the 
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V 


maximum  lies  closer  to  the  axis  [cf.  Figures  (III. A. 2)  and  (III. A. 3)]. 
Hence,  the  more  the  pressure  falls  in  the  core,  the  more  fluid  sinks  into 
the  boundary  layer  to  spiral  inward,  erupt  upward,  and  cause  further 
pressure  reduction  in  the  core.  If  the  crucial  early  competition  is 
resolved  in  favor  of  the  organized  convection,  ultimately  the  particles 
erupting  out  of  the  boundary  layer  rise  so  quickly  that  they  lie  on  a 
moist  adiabat,  and  the  greatly  lightened  core  is  entirely  flushed  of  its 
original  fluid. 

No  mention  has  yet  been  made  of  the  eye.  This  is  in  direct  contrast 
with  the  description  of  intensification  given  in  Palmen  and  Newton  (1969) 
in  which  the  eye  is  depicted  as  being  gradually  formed  as  the  pressure 
deficit  develops.  In  the  Carrier  model  the  central  core  is  completely 
flushed  of  ambient-like  air,  so  that  the  air  in  the  core  lies  on  a  moist 
adiabat  based  on  sea-level  ambient  conditions,  before  any  trace  of  an  eye 
is  to  be  found.  A  Rankine-vortex-like  swirl  holds  everywhere.  From  this 
fully  developed  one-cell  structure,  a  two-cell  structure  with  a  calm 
center  region  emerges  rapidly,  probably  in  much  less  than  an  hour,  owing 
to  inertial  oscillation,  in  the  following  way. 

As  the  pressure  falls  in  the  core  relative  to  ambient,  the  particles 
in  I  necessarily  move  in  closer  to  the  axis  to  permit  a  compensating 
centrifugal  force  to  develop.  Once  the  core  is  flushed  and  moist  adiabatic 
ascent  characterizes  the  full  height  of  the  core,  no  further  pressure 
deficit  can  be  generated.  By  inertia,  the  spinning  particles  continue  to 
move  in,  a  dynamic  imbalance  is  created,  and  a  radial  acceleration  develops 
to  force  the  particles  away  from  the  axis  of  symmetry.  This  reverse  motion 
creates  a  rarefaction  at  the  center,  and  relatively  dry  warm  motionless  air 
sinks  down  the  axial  column  to  form  an  eye.  This  air  may  be  air  from 
above  the  storm  or  rained-out,  slowly  swirling  air  entrained  out  of  the 
top  of  the  moist  adiabatic  column  [cf.  Figure  (III. A. 4)].  Because  there 
is  no  appreciable  swirl  (hence  no  associated  pressure  gradient)  in  the  eye, 
there  is  no  frictional  boundary  layer  under  the  eye.  The  moist  adiabatic 
column  becomes  an  annulus  displaced  from  the  axis,  i.e.,  the  eyewall;  the 
inertial  oscillations  of  the  eyewall  eventually  damp  in  time. 
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Fig.  III. A. 4  Schematic  picture  of  the  flow  configuration  which 
prevails  when  R  is  increasing  and  a  nascent  eye  is  being  filled  with 
relatively  dry  and  motionless  air,  which  sinks  down  from  the  top  of  the 
storm  under  dry-adiabati c  recompression.  With  the  formation  of  an  eye 
lighter  in  weight  than  the  eyewall,  the  terminal  stages  of  intensification 
and  the  beginnings  of  quasi  steady  mature-stage  structure  are  realized. 


